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ABSTRACT 



Type I X-ray bursts are violent stellar events that take place on the H/He- 
rich envelopes of accreting neutron stars. We have investigated the role played 
by uncertainties in nuclear processes on the nucleosynthesis accompanying these 
explosive phenomena. Two different approaches have been adopted, in the frame- 
work of post-processing calculations. In the first one, nuclear rates are varied 
individually within uncertainties. Ten different models, covering the characteris- 
tic parameter space for these stellar events, have been considered. The second, 
somewhat complementary approach involves a Monte Carlo code in which all 
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nuclear rates are randomly varied within uncertainty limits simultaneously. All 
in all, about 50,000 post-processing calculations, with a network containing 606 
nuclides (H to 113 Xe) and more than 3500 nuclear processes, have been performed 
in this work. A brief comparison between both procedures is outlined together 
with an overall account of the key nuclear reactions whose uncertainties have the 
largest impact in our X-ray burst nucleosynthesis studies. 

Subject headings: Nuclear reactions - nucleosynthesis - abundances - stars: neu- 
tron - X-rays: bursts 



1. Introduction 

Type I X-ray bursts (hereafter, XRBs) are cataclysmic stellar events powered by ther- 
monuclear runaways (TNR) in the H/He-rich envelopes accreted onto neutron stars in close 
binary systems (see reviews by Bildsten 1998; Lewin et al. 1993, 1995; Psaltis 2006; Schatz 
& Rehm 2006; Strohmayer & Bildsten 2006). They constitute the most frequent type of 
thermonuclear stellar explosion in the Galaxy (the third, in terms of total energy output 
after supernovae and classical novae) because of their short recurrence period (hours to 
days). More than 80 Galactic low-mass X-ray binaries exhibiting such bursting behavior 
(with r burst ~ 10 - 100 sec) have been found since the discovery of XRBs by Grindlay et al. 
(1976), and independently, by Belian et al. (1976). 

Modeling of type I X-ray bursts and their associated nucleosynthesis (see pioneering 
models by Woosley & Taam 1976, Maraschi & Cavaliere 1977, and Joss 1977) has been 
extensively addressed by different groups, reflecting the astrophysical interest in determining 
the nuclear processes that power the explosion as well as in providing reliable estimates for 
the chemical composition of the neutron star surface. Indeed, several thermal (Miralda- 
Escude, Paczyhski, & Haensel 1990; Schatz et al. 1999), electrical (Brown & Bildsten 1998; 
Schatz et al. 1999), and mechanical properties (Bildsten & Cutler 1995; Bildsten & Cumming 
1998) of the neutron star depend critically on the specific chemical abundance pattern of 
its outer layers. The diversity of shapes in XRB light curves (L peak ~ 3 x 10 38 erg s _1 
- see e.g. Galloway et al. (2007), Lewin et al. (1993), and Kuulkers et al. (2003)) is 
also likely due to different nuclear histories (see Heger et al. 2008, for an account of the 
interplay between long bursts and the extension of the rp-process in XRBs), suggesting that 
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the final composition is not unique. It is worth noting that, as discussed by Taam (1980) 
and Woosley et al. (2004), the properties of the bursts recurring in a given stellar source are 
affected by compositional inertia; that is, they are sensitive to the fact that accretion proceeds 
onto the ashes of previous bursts. Indeed, this compositional inertia seems to reduce the 
expected recurrence times between bursts (particularly for scenarios involving accretion of 
low-metallicity matter). Moreover, these ashes may provide characteristic signatures such as 
gravitationally redshifted atomic absorption lines from the surface of the neutron star, which 
could be identified through high-resolution X-ray spectra (see Cottam, Paerels, & Mendez 
2002; Bildsten, Chang, & Paerels 2003; Chang, Bildsten, & Wasserman 2005; Chang et 
al. 2006; Weinberg, Bildsten, & Schatz 2006), providing a valuable tool to constrain XRB 
models. Furthermore, the mechanism responsible for superbursts, highly energetic (~ 10 42 
erg), long duration (r ~ 10 4 s) bursts (Cornelisse et al. 2000; Wijnands 2001, and references 
therein), first reported from the source 4U 1735-44 (Cornelisse et al. 2000), also depends upon 
XRB nucleosynthesis. These phenomena are attributed to ignition in a C-layer accumulated 
from successive type I XRBs (first proposed by Woosley & Taam 1976; see also, Taam & 
Picklum 1978, Brown & Bildsten 1998, Cumming & Bildsten 2001, Schatz, Bildsten, & 
Cumming 2003, Weinberg, Bildsten, & Brown 2006, and Weinberg & Bildsten 2007). The 
potential impact of XRB nucleosynthesis on Galactic abundances is still a matter of debate: 
although ejection from a neutron star is unlikely because of its large gravitational potential 
(matter accreted onto a neutron star of mass M and radius R releases GMm p /R ~ 200 
MeV/nucleon, whereas only a few MeV/nucleon are released from thermonuclear fusion), 
radiation-driven winds during photospheric radius expansion may lead to ejection of a tiny 
fraction of the envelope (containing nuclear processed material -see Weinberg et al. 2006a 
and MacAlpine et al. 2007). Indeed, XRBs may help to explain the Galactic abundances of 
the problematic light p-nuclei (Schatz et al. 1998). 

With a neutron star as the underlying compact object, temperatures and densities in 
the accreted envelope reach quite high values: T peak > 10 9 K, and p ~ 10 6 g.cm -3 (note that 
during superbursts, however, densities may exceed 10 9 g.cm -3 -see Cumming & Bildsten 
2001). These values are at least an order of magnitude larger than in a typical classical nova 
outburst - see e.g. Jose & Hernanz (1998). As a result, detailed nucleosynthesis studies 
require the use of hundreds of isotopes, up to the SnSbTe mass region (Schatz et al. 2001) 
or beyond (the nuclear activity in the XRB nucleosynthesis studies of Koike et al. (2004) 
reaches 126 Xe), and thousands of nuclear interactions. The main nuclear reaction flow is 
driven by the rp-process (rapid proton-captures and /5 + -decays), the 3a-reaction and the 
ap-process (a sequence of (a,p) and (p,7) reactions), and is expected to proceed far away 
from the valley of stability, merging with the proton drip-line beyond A = 38 (Schatz et al. 
1999). 
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Until recently, because of computational limitations, studies of XRB nucleosynthesis 
have been performed using limited nuclear reaction networks, truncated around Ni (Woosley 
& Weaver 1984; Taam et al. 1993; Taam, Woosley, & Lamb 1996 -all using a 19-isotope 
network), Kr (Hanawa, Sugimoto, & Hashimoto 1983 -274-isotope network; Koike et al. 
1999 -463 nuclides), Cd (Wallace & Woosley 1984 -16-isotope network), or Y (Wallace & 
Woosley 1981 -250-isotope network). On the other hand, Schatz et al. (1999, 2001) have 
carried out very detailed nucleosynthesis calculations with a network containing more than 
600 isotopes (up to Xe, in Schatz et al. 2001), but using a one-zone approach. Koike et al. 
(2004) have also performed detailed one- zone nucleosynthesis calculations, with temperature 
and density profiles obtained from a spherically symmetric evolutionary code, linked to a 
1270-isotope network extending up to 198 Bi. 

Note however that different numerical approaches and approximations (hydrodynamic 
simulations with limited networks or one- zone calculations with detailed networks) have been 
adopted in all those works, and hence, the predicted nucleosynthesis in each case has to be 
taken with caution. Indeed, recent attempts to couple hydrodynamic stellar calculations 
(in 1-D) and detailed networks include Fisker et al. (2004, 2006, 2007, 2008) and Tan et 
al. (2007) (using networks of ~ 300 isotopes, up to 107 Te), Jose & Moreno (2006) (using 
a network of 2640 nuclear reactions and 478 isotopes, up to Te) and Woosley et al. (2004) 
(using up to 1300 isotopes with an adaptive network). 

Most of the nuclear reaction rates required for these extensive nucleosynthesis calcu- 
lations (in particular, those related to heavy species) rely on theoretical estimates from 
statistical models, and therefore may be affected by significant uncertainties. To date, only 
partial efforts have been made to quantify the impact of such nuclear uncertainties on the 
overall abundance pattern accompanying XRBs (Wallace & Woosley 1981; Schatz et al. 
1998; Iliadis et al. 1999; Koike et al. 1999, 2004; Thielemann et al. 2001; Fisker et al. 2004, 
2006, 2008; Amthor et al. 2006), revealing a complex interplay between nuclear activity 
and the shape of the light curve (Hanawa et al. 1983, Woosley et al. 2004). It has to be 
stressed, however, that simple, parametrized models cannot adequately address this issue, 
since nuclear changes in the innermost envelope shell do not show up immediately at the 
surface of the star. 

But in order to provide reliable estimates of the composition of neutron star envelopes 
it is of paramount importance to fully identify the key reactions whose uncertainties have 
the largest impact on XRB yields (reactions that probably deserve further improvement 
through nuclear physics experiments at dedicated facilities). To achieve such a goal, we have 
performed a comprehensive study of the effects of thermonuclear reaction-rate variations on 
type I X-ray burst nucleosynthesis, sampling the overall parameter space. The scale of this 
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study makes state-of-the-art hydro dynamic models computationally prohibitive. Hence, we 
rely on post-processing calculations, coupling a detailed and fully updated nuclear reaction 
network with temperature and density profiles extracted from literature. 

The structure of the paper is as follows. In Section 2, we outline the two approaches 
adopted in this paper, in the framework of post-processing calculations with temperature 
and density profiles extracted from the literature. In Section 3, we extensively describe the 
results obtained from the different models, through a sensitivity study in which all reaction 
rates are varied individually within uncertainty limits. This is compared with Monte Carlo 
simulations in which all nuclear reactions are simultaneously varied with randomly sampled 
enhancement factors (Section 4). A summary of the main conclusions reached in this study, 
including a list of key reactions whose uncertainties deeply influence the final type I XRB 
yields, is presented in Section 5. 

2. Models & Input Physics 

2.1. Rate uncertainties and post-processing calculations 

Nucleosynthesis studies in XRB conditions are computationally challenging because of 
the characteristic timescales and peak temperatures achieved, the large number of nuclear 
species involved (several hundred), and the huge networks (with several thousand nuclear 
processes) required to address the rp- and ap-processes occuring for H- and He-rich mixtures 
accreted onto the surface of a neutron star. 

Two different, somewhat complementary approaches, based on post-processing calcula- 
tions with temperature and density profiles, can be adopted in a comprehensive sensitivity 
study. In the first one, all rates are varied individually (see Section 3) within uncertainty 
limits so as to check the impact of each nuclear process on the final yields. At least 2 post- 
processing calculations are required per nuclear rate to account for the upper and lower limits 
posed by its associated uncertainty. For a network containing several thousand interactions 
the overall number of post-processing calculations is indeed extraordinarily large. It is worth 
mentioning that this technique has been previously applied to a large number of astrophysi- 
cal sites, including nucleosynthesis in the Sun (Bahcall et al. 1982), type II supernovae (The 
et al. 1998; Jordan, Gupta, & Meyer 2003), classical nova explosions (Iliadis et al. 2002), 
primordial (Big Bang) nucleosynthesis (Coc et al. 2002, 2004), intermediate-mass AGB stars 
(Izzard et al. 2007), and Type I X-ray bursts (Amthor et al. 2006). 

A second approach is based on Monte Carlo techniques (see Section 4). Here, random 
enhancement factors (often adopted to follow a log-normal distribution) are applied to each 
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nuclear process of the network simultaneously. The impact on the final yields is then tested 
through a series of post-processing calculations. This approach requires a large number of 
trial simulations in order to be statistically sound and has been already applied to Big Bang 
nucleosynthesis studies (Krauss & Romanelli 1990; Smith et al. 1993), nova nucleosynthesis 
(Smith et al. 2002; Hix et al. 2002, 2003) and also to Type I X-ray bursts (Roberts et al. 
2006). 

Both approaches have been adopted in this study (although more emphasis has been 
placed on the first one). An interesting issue, recently raised by Roberts et al. (2006), is the 
feasibility of the first method, as compared with the Monte Carlo approach, to properly ad- 
dress the higher-order correlations between input rates and XRB model predictions because 
of the large number of reactions simultaneously involved in the production and destruction 
of each element. It is also the goal of this manuscript to test this conjecture by comparing 
the results from an individual- variation sensitivity study with those obtained with a Monte 
Carlo approach. 



2.2. Initial Models: Temperature and Density profiles 

In this work, we have used three temperature- density versus time profiles directly ex- 
tracted from the literature (see Fig. 1). Model K04 (Koike et al. 2004) is based on a 
spherically symmetric model of accretion onto a 1.3 M neutron star (with Rns=8.1 km). 
This model is characterized by a peak temperature of 1.36 GK, densities ranging between 
(0.54 - 1.44) xlO 6 g.cm~ 3 , and a burst duration of ~ 100 s. The initial composition is taken 
directly from Koike et al. (2004), with X=0.73, Y=0.25 and Z=0.02 -roughly solar. 

Model F08 (Fisker et al. 2008) is based on 1-D general relativistic hydrodynamic cal- 
culations. It achieves a peak temperature of only 0.993 GK, densities ranging between (2.07 
- 5.14) xlO 5 g.cm -3 , and a burst duration of ~ 50 s. The initial composition is taken di- 
rectly from the hottest burning zone of the hydrodynamic models of Fisker et al. (2008), 
with X=0.40, Y=0.41 and Z=0.19 -a metallicity about ten times solar; see the original 
manuscript for a detailed account of the distribution of those metals. 

Finally, Model S01 (based on Schatz et al. 2001) achieves the largest temperature of all 
models extracted from the literature, with T pea k = 1.907 GK. The shape of this temperature 
profile, based on one-zone model calculations, is quite different from those corresponding to 
Models K04 and F08 (which approximately follow an exponential decay): indeed, it shows a 
long plateau at about T ~1.4 GK before the final decline. Since no density profile is directly 
available for this model, and considering the relatively small variation in density during an 
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XRB, we have scaled K04 to match the values reported in Schatz et al. (2001), resulting in 
densities ranging between (0.54 - 1.73) xlO 6 g.cm" 3 . The duration is ~ 300 s. The initial 
composition is X=0.718, Y=0.281 and Z=0.001 -a metallicity about 20 times lower than 
solar. Because of the lack of information on the specific distribution of the metallicity, we 
assume that all metals correspond to 14 N (see also Woosley et al. 2004), following the rapid 
rearrangement of CNO isotopes that naturally occurs early in the burst. 

The three models mentioned above partially cover the parameter space in XRB nucle- 
osynthesis calculations. However, it is difficult, if not impossible, to disentangle the specific 
role played by the initial metallicity, the peak temperature achieved and the duration of 
the bursts, since all these models are characterized by different values for these physical 
quantities. Hence, to appropriately address this issue, we have generated additional mod- 
els through parameterization of the Koike et al. (2004) temperature-density-time profiles. 
Parameterized models have been used previously to study, for example, nucleosynthesis in 
novae (Boffin et al. 1993), the alpha-rich freeze-out in Type II supernovae (The et al. 1998, 
Jordan et al. 2003) and XRBs (Wallace and Woosley 1981, 1984). 

To evaluate the role played by the duration of the burst (taken as the characteristic 
timescale of the temperature and density profiles) in the extension of the nuclear path, we 
have scaled the K04 profiles in duration by a factor 0.1 (short burst: Model K04-B1) and 10 
(long burst: Model K04-B2) while preserving Tpeak and the initial chemical composition. 

The role played by the initial metallicity (reflecting that of the stream of infalling ma- 
terial during accretion) has been tested through 3 models: a low- metallicity model (Z=10~ 4 ; 
Model K04-B3), and two high metallicity ones (with Z=0.19; Models K04-B4 & K04-B5). 
The initial composition of Models K04-B3 and K04-B4 is scaled from the distribution re- 
ported in Koike et al. (2004). In contrast, for Model K04-B5, we have adopted the distribu- 
tion given in Fisker et al. (2008). This will also allow us to test if a different distribution of 
metals has an impact on the final yields. 

Finally, two additional models have been constructed to test the effect of varying the 
peak temperature achieved during the explosion: we have scaled the temperature and density 
versus time profiles from Koike et al. (2004) to attain peak values of 0.9 x 10 9 K (lower peak 
temperature; Model K04-B6) and 2.5 x 10 9 K (higher peak temperature; Model K04-B7), 
while preserving the initial composition and the duration of the burst. 

It is worth noting that post-processing calculations are not suited to derive absolute 
abundances (or, as mentioned above, to provide any insight into light-curve variations and 
energetics) since they rely only on temperature and density versus time profiles evaluated 
at a given location of the star (usually, the innermost shells of the envelope). Indeed, it is 
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likely that the evolution at other depths will be characterized by a different set of physical 
conditions. Furthermore, adjacent shells will eventually mix when convection sets in, altering 
the chemical abundance pattern in those layers. However, this approach is reliable to identify 
the key processes governing the main nuclear activity at the specific temperature and density 
regimes that characterize such bursting episodes. Hence, the goal of this paper is to provide 
the reader with a list of key nuclear processes whose uncertainties have the largest influence 
on the final yields, covering as much as possible the proper range of temperatures, densities, 
and timescales that characterize XRBs. Any attempt to properly quantify the extent of this 
impact (for instance, in terms of absolute abundances) however, must rely on state-of-the-art 
hydrodynamic codes coupled to detailed nuclear networks. 



2.3. Nuclear Reaction Network 

The temperatures and densities achieved in XRB nucleosynthesis are sufficiently high 
so that many nuclear reactions, especially those with relatively small Q-values, achieve an 
equilibrium between the forward and reverse process. For example, if the reaction A(p, j)B 
has a small Q-value, then the strong B(j,p)A photodisintegration will give rise to a small 
equilibrium abundance of nuclei B that may then capture another proton. Such cases, 
which are called (sequential) two-proton captures, must be considered carefully since they 
represent waiting points (and may even be candidates for termination points) for a continuous 
abundance flow toward heavier-mass nuclei. Among the most important of these waiting 
points are 64 Ge, 68 Se and 72 Kr. When a reaction rate equilibrium has been established, the 
most important nuclear physics information needed is the reaction Q-value (which enters 
exponentially in the effective decay constant of A) of the link A(p, ^)B and the reaction 
rate of the second-step reaction B(p, 7)C, but the rate of the reaction A(p, ^)B is irrelevant. 
Note that for the three nuclei mentioned above, the (p, 7) Q-values have not been measured 
directly yet. According to Audi et al. (2003a), their predicted values amount to Q — 
—80 ± 300 keV, —450 ± 100 keV and —600 ± 150 keV, respectively. Some encouraging 
progress has been made through mass measurements of 64 Ge (Clark et al. 2007, Schury 
et al. 2007), 68 Se (Clark et al. 2004, Wohr et al. 2004, Chartier et al. 2005), and 72 Kr 
(Rodriguez et al. 2004). 

In order to set up the reaction library for our study, we started by adopting the proton 
drip line from Audi et al. (2003a,b). For reactions A(p, 7)5 with Q-values below 1 MeV the 
rates (below the element Pd) are calculated using the Hauser-Feshbach code MOST (Goriely 
1998; Arnould & Goriely 2006) and the corresponding reverse photodisintegrations are com- 
puted with the extrapolated proton separation energies from Audi et al. (2003a). Notice that 
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nuclear masses enter twice: first in the calculation of the forward Hauser-Feshbach rate, and 
second in the calculation of the reverse photodisintegration rate. The second-step reactions, 
B(p, 7)C, are also computed with the code MOST. We had to include many nuclides beyond 
the proton drip line in order to properly account for sequential two-proton captures. For all 
other reactions for which experimental rates are not available, we used the results from the 
Hauser-Feshbach code NON-SMOKER (Rauscher & Thielemann 2000). 

Our entire nuclear network consists of 606 nuclides ranging from 1 H to 113 Xe. All 
charged-particle induced reactions between these isotopes, and their corresponding reverse 
processes, have been included. Experimental information was only available for the rates of 
a small subset of reactions. These rates are adopted from Angulo et al. (1999), Iliadis et al. 
(2001), and some recent updates for selected reactions. Neutron captures are disregarded 
since our early test calculations have revealed that they play a minor role in XRB nucleosyn- 
thesis. All reaction rates incorporate the effects of thermal excitations in the target nuclei 
(Rauscher & Thielemann 2000). 

For the weak interactions, we use laboratory decay rates (Audi et al. 2003b); the impact 
of /3-delayed nucleon emission has also been considered. For a discussion of employing stellar 
versus laboratory decay rates, see Woosley et al. (2004). Also, note that many computed 
stellar decay rates (Fuller, Fowler & Newman 1982a,b) do not converge to their laboratory 
values at lower temperatures and densities. In order to avoid this problem we decided 
against using stellar decay rates. Although the use of properly converging stellar decay rates 
is in principle preferred, we doubt that the effect would be large in XRB nucleosynthesis 
calculations. 



3. Individual Reaction-Rate Variations 

Some general procedures must be outlined before a detailed analysis of the main results is 
presented. For all 10 models considered in this Section, a first calculation, with standard rates 
(as described in Section 2.3), has been performed. This is used to scale the level of changes 
(indicated in Tables 1-10, & 13-14 as Xi/X itSt d, that is, the ratio of mass fractions obtained 
with a modified network -resulting from our exploration of uncertainties- to those obtained 
with our standard network). Following this, a whole series of post-processing calculations 
have been computed, in which each nuclear rate is varied individually by a factor 10 up and 
down, for each of the 10 models. 

Several important remarks have to be made at this stage: first, we have chosen a factor 
of 10 for the level of uncertainty affecting theoretical reaction rate estimates, in general. 
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Other authors (see Schatz 2006, and Amthor et al. 2006) claim, instead, that excitation 
energies of theoretically calculated levels for XRB conditions may suffer uncertainties of ~ 
100 keV, which translates into an overall uncertainty in some rates that may reach several 
orders of magnitude. It is also worth noting that efforts to ascertain systematic uncertainties 
through a direct comparison between different theoretical models (see Section 5.1) often 
report variations smaller than a factor of 10. Second, some nuclear rates, such as the triple- 
alpha reaction, some (p, 7) reaction rates on low-mass targets, or most of the /3-decay rates 
(note again that we used ground-state or laboratory /3-decay rates), are known with much 
better precision. This has been addressed when necessary and is discussed in the following 
sections of this manuscript. Third, it is important to stress that forward and reverse reactions 
have been varied simultaneously with the same uncertainty factor, to preserve the principle 
of detailed balance. Fourth, a particular strategy has been adopted to evaluate the impact 
of /3-decay rate uncertainties: they have initially been varied by a factor of 10, up and down 
as a way to discriminate the key /5-decay rates from the rest (see Table 11). Then, additional 
post-processing calculations focused on these key rates have been performed assuming more 
realistic uncertainties (around 30%, according to ground-state /3-decay rate uncertainties; see 
Audi et al. 2003b). Fifth, some additional tests aimed at identifying the influence of Q- value 
variations have been performed using the uncertainties estimated by Audi et al. (2003a) (see 
Tables 12-14). This part of our study has been restricted to Models K04 and F08, to 
proton-capture reactions on A < 80 nuclei, and to reactions with \Q\ < 1 MeV, for which 
the estimated uncertainty exceeds 50 keV (see Table 12). Sixth, since the main goal of this 
paper is to identify the key nuclear processes whose uncertainties have the largest impact 
on XRB nucleosynthesis, we have restricted the analysis (in the Tables and forthcoming 
discussion) to nuclear species which achieve a mass fraction of at least 10" 5 at the end of 
the burst, and deviate from the abundances computed with standard rates by at least a 
factor of 2. And seventh, for the sake of clarity and conciseness, the isotopes displayed and 
discussed throughout this manuscript correspond to species that are either stable or have a 
half-life longer than 1 hour (the rest are assumed to fully decay at the end of the burst and 
consequently, are added to the corresponding stable or long-lived daugther nuclei). 



3.1. Results 



In this Section, we report results from a series of ~40,000 post-processing calculations 
(requiring 14 CPU-months) performed for this study of X-ray burst nucleosynthesis. 
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3.1.1. Model K 04 

Here, we will describe in detail the results obtained from our post-processing calculations 
using temperature and density versus time profiles extracted from Koike et al. (2004) [Model 
K04]. A general discussion of the main nuclear path achieved for these profiles can be found 
in Iliadis (2007). Our main results, from individually varying each reaction rate by a factor 
10 (up and down), are summarized in Table 1. 

Because of the moderate peak temperature achieved, the extent of the nuclear activity 
(defined as the heaviest isotope with Xj > 10~ 2 , in the final yields) reaches 96 Ru. The most 
abundant species (stable or with a half-life > 1 hr) at the end of the burst are H (0.20, by 
mass), 4 He (0.021), 68 Ge (0.20), 72 Se (0.13), 64 Zn (0.071), and 76 Kr (0.074). Qualitatively, 
this is in agreement with the results reported from Koike et al. (2004) (Model 2, Table 8), in 
which the most abundant species, at the end of the burst, are also 68 Ge (0.17), 72 Se (0.14), 
64 Zn (0.27), and Kr (0.078). Note that the final mass fraction of 4 He reported by Koike et 
al. (2004), 0.011, is also comparable to the value reported in this paper. However, Koike's 
model is fully depleted of H (9.76 x 10~ 17 , by mass), whereas some H remains at the end of 
our simulations. This could result from differences in the adopted nuclear reaction networks 
as well as to the hybrid use of evolutionary and post-processing calculations in Koike et al. 
(2004). This discrepancy is of some importance as the presence of H in the ashes may induce 
marginal nuclear activity driving a second burst if the set of necessary conditions for a TNR 
are satisfied (see Woosley & Weaver 1984, for details). 

Our study reveals that the most influential reaction, by far, is the 3 a, which affects a 
large number of species ( 4 He, 18 F, 21 Ne, 24 > 25 Mg, 28 " 30 Si, 33 < 34 S, 36 " 38 Ar, 41 Ca, ^^Ti, 49 V, 
50 Cr, 53 Mn, 54 Fe, 57 " 59 Ni, 61 ' 63 Cu, 62 Zn, 102 Pd, and 103 . 104 Ag), when its nominal rate is varied 
by a factor of 10, up and down. A similar result on the importance of the 3 a reaction has 
been previously reported by Roberts et al. (2006), using a Monte Carlo approach. But, as 
discussed above, this uncertainty factor is far too large and accordingly, this reaction has 
been removed from Tables 1-10. Hereafter, we will drop from the discussion of the different 
models any additional comment on the importance of the 3a reaction which, when arbitrarily 
varied by a factor of 10, systematically becomes the single most influential reaction of the 
whole network, for all 10 models. Indeed, additional tests performed with a more realistic 
uncertainty of ±40% (see Angulo et al. 1999) show no effect on any individual isotope, for 
any of the 10 models (the impact of nuclear uncertainties affecting the triple alpha rate on 
the total energy output will be specifically addressed in Section 5). 



2 Note that according to Tur et al. (2006) the 3a rate is known to ± 12%. 
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From the several thousand nuclear processes, we find only 56 reactions-and the corre- 
sponding reverse reactions- have an impact on the final yields when their rates are varied 
by a factor of 10, up and down for this model (see Table 1). Furthermore, our study reveals 
that the impact of most of these reactions is restricted to the vicinity of the target nuclei. 
A clear example is 15 0(a,7): when its nominal rate is multiplied by a factor 0.1, only one 
isotope, 15 N, is modified by more than a factor of 2 (as compared with the mass fraction 
obtained with the recommended rate -but see also Section 5). Indeed, the most influential 
reaction is by far 65 As(p,7), perhaps expected due to its bridging effect on the 64 Ge- waiting 
point. To a lesser extent, 96 Ag(p,7), and 102 In(p,7) also show an impact on a number of 
nuclear species. 

The most important /3-decay rates identified in this Model are those of 64 Ge, 68 Se, 
72 Kr, 76 Sr, 80 Zr, 88 Ru, 92 Pd, and "In (see Table 11). But it is important to stress that 
additional calculations, in which these (laboratory) /3-decay rates were varied within realistic 
uncertainties (half-lives varied by ~ ±30%), have revealed no effect on any single isotope 
(see Woosley et al. 2004, for a sensitivity study based on variations of groups of weak rates, 
including all positron emission rates for nuclei heavier than 56 Ni, by an order of magnitude). 

From the list of nuclear reactions selected for the Q-value variation study (see Table 

12, and the discussion in Section 3), only 26 P(p, 7) 27 S, 46 Cr(p, 7) 47 Mn, 55 Ni(p, 7) 56 Cu, 
60 Zn(p, 7) 61 Ga, and 64 Ge(p, 7) 65 As show some impact on the final yields when their De- 
values are varied between one sigma uncertainty bounds. However, as summarized in Table 

13, the effects are restricted to single isotopes, except for the last reaction, 64 Ge(p, 7), whose 
influence ranges between 64 Zn and 104 Ag. Indeed, the significance of this reaction on the 
64 Ge waiting point is well-known (see e.g., Woosley et al. 2004, Thielemann et al. 2001, 
Schatz 2006, Fisker et al. 2008, and references therein). Mass measurements on 64 Ge have 
indeed been performed (Clark et al. 2007; Schury et al. 2007). When combining these 
measurements with the 65 As mass value given by Audi et al. (2003a), the corresponding one 
sigma uncertainty limits in the 64 Ge(p, 7) 65 As Q-value (—74 ±303 keV, and —46 ±302 keV, 
respectively) are well covered by the range used in our study (Q = —80 ± 300 keV; Audi 
et al. (2003a)). Clearly, mass measurements on 65 As are still required to better determine 
this Q-value. These are challenging experiments because of the difficulty of producing 65 As 
(Clark, private communication). 



3.1.2. Model S01 



Here, we summarize the main results obtained from our post-processing calculations 
with the temperature versus time profile extracted from the one- zone model of Schatz et al. 
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(2001) [Model SOI]. 

Because of the larger peak temperature achieved in this model, the main nuclear activity 
extends to heavier nuclei ( 107 Cd, for X > 10~ 2 ) than in Model K04. The most abundant 
species (stable or with a half-life > 1 hr) at the end of the burst are now H (0.071, by mass), 
4 He (0.013), and the heavy isotopes 104 Ag (0.328), 106 Cd (0.244), 103 Ag (0.078), and 105 Ag 
(0.085). This chemical pattern bears some resemblance to that reported by Schatz et al. 
(2001), and also to a similar model (zM) computed by Woosley et al. (2004) with a 1-D 
hydrodynamic code. For example, a large overproduction of the p-nucleus 106 Cd is explicitly 
mentioned in Schatz et al. (2001), Fig. 4; as well, Woosley et al. (2004) (Fig. 7) find a large 
mass fraction of 106 Sn (which will decay into 106 Cd), followed by 64 Ge, 68 Se, 104 In (which 
will decay into 104 Ag), 4 He, and H. It is worth noting that the final amount of H is much 
smaller than in Model K04: here, H and He are almost depleted and hence the next burst 
will necessarily require the piling-up of fresh H/He-rich fuel on top of this H/He-depleted 
shell. 

Only 64 reactions (and the corresponding reverse processes) from the overall nuclear 
network turn out to have an impact on the final yields when their rates are varied by a 
factor of 10, up and down (Table 2). As reported for Model K04, the influence of these 
reactions is often limited to the vicinity of the target nuclei. Remarkable exceptions are 
56 Ni(o;,p), 59 Cu(p,7), and to a lesser extent, 103 Sn(o;, p). Notice also that, among the most 
important reactions, only three a-captures ( 15 0(a,7), 56 Ni(o;,p), and 103 Sn(o;, p)) are found. 
The role played by 15 0(a,7) is as marginal as for Model K04 however (but see discussion in 
Section 5). 

Concerning the most influential /3-decay rates, we found similarities with the results 
reported for Model K04. Hence, the /9-decay rates of 68 Se, 72 Kr, 80 Zr, 88 Ru, and 92 Pd, are 
important, but due to the higher peak temperatures reached in this Model, the list also 
extends to heavier species, including 101 ~ 106 Sn and 106 Sb (see Table 11). Additional studies 
of the specific impact of these decays with realistic uncertainty bounds (Audi et al. 2003b) 
have also been performed for this case, revealing no effect on any particular isotope. 



3.1.3. Model F0 8 

Since this model achieves the lowest peak temperature of the three models extracted 
from the literature, its main nuclear path is somewhat more limited, reaching only (X > 10~ 2 ) 
72 Se: the most abundant species (stable or with a half-life > 1 hr) at the end of the burst 
are now 4 He (0.085), 56 Ni (0.13), 60 Ni (0.38), 12 C (0.040), 28 Si (0.041), and 64 Zn (0.034). 
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Moreover, H has been depleted down to 10~ n by mass. Qualitatively, there is good agreement 
with the relevant chemical pattern reported in Fisker et al. (2008), Fig. 20. It is worth noting 
that the large amount of 12 C obtained at the end of the burst may have implications for the 
energy source that powers superbursts (see Sect. 1). 

Again, only a few reactions -62 and the corresponding reverse reactions- have an impact 
on the final yields when their rates are varied by a factor of 10, up and down (Table 3). The 
most influential reactions are 12 C(a,7), 269 Al(o;, p), 57 Cu(p,7), 61 Ga(p,7), and 85,86 Mo(p,7). 
Additional discussion and tests on the impact of the 12 C(a,7) rate are presented in Section 
5. Thrirteen a/p/ia-capture reactions show some impact on the yields, a larger number than 
in K04 and SOI. Indeed, all key reactions affecting A < 30 species are mainly (a,j) or (a, p) 
reactions. It is worth mentioning, however, that 15 0(a,7) is absent from this list (see Sect. 
5). 

The most influential /5-decay rates are now shifted towards lighter unstable isotopes, 
because of the lower T peak value achieved. These include decays of 26 ' 29 Si, 29 ' 30 S, 33 ' 34 Ar, 
37 > 38 Ca, 60 Zn, 68 Se, 72 Kr, 76 Sr, and 80 Zr (Table 11). As in K04 and S01 though, no effect is 
seen on any isotope when realistic uncertainties are adopted for these decay rates. 

Finally, the only nuclear reaction from Table 12 that has a significant impact on the 
final yields when its Q-value is varied between one sigma uncertainty bounds is 45 Cr(p, 7) 
(see Table 14). 



3.1.4. Effect of the duration of the burst: Model K04-B1 vs. K04-B2 

Because of the shorter duration of Model K04-B1, as compared with Model K04, the 
main nuclear path extends up to 68 Ge. In contrast, the longer burst of Model K04-B2 drives 
the main path up to 106 Cd. 

The most abundant species at the end of the bursts are: H (0.42, by mass), 4 He (0.061), 
60 Ni (0.10), 64 Zn (0.31), and 68 Ge (0.023), for the shorter burst of Model K04-B1; and 4 He 
(0.0051), 68 Ge (0.159), 72 Se (0.121), 76 Kr (0.077), 80 Sr (0.040), 90 Mo (0.042), 94 Mo (0.055), 
98 Ru (0.040), 102 Pd (0.034), 103 Ag (0.060), 104 Ag (0.097), and 105 Ag (0.042), for Model K04- 
B2. The final amount of H in the latter model has decreased down to 10 -15 . The larger 
extension of the main nuclear activity reported for Model K04-B2 is a direct consequence 
of the longer exposure times to higher temperatures. However, the depletion of H results 
from the complex interplay between exposure time to high temperatures, which decreases 
the overall H content, and the effects of photodisintegrations, which raise the H abundance. 
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Concerning nuclear uncertainties, Model K04-B1 is characterized by 28 critical reactions, 
the most important ones being 61 Ga(p,7) and 65 As(p,7), with a marginal role played by 
18 Ne(o;,p) and 31 Cl(p,7) (Table 4). For Model K04-B2, with 51 critical reactions, the most 
important are 65 As(p,7) and 32 S(o;,7), followed, to some extent, by 12 C(a, 7), 61 Ga(p, 7), 
75 Rb(p, 7), 84 Zr(p, 7), 92 Ru(p, 7), 93 Rh(p, 7), and 96 Ag(p, 7) (Table 5). Only 2 a — capture 
reactions appear to be influential in Model K04-B1 (one being 15 0(a,7), see Table 4), whereas 
uncertainties affecting 12 a-capture reactions turn out to be critical for Model K04-B2 (see 
Table 5). 

Table 11 lists the most important /3-decay rates for these models, when these rates are 
varied by a factor of 10 up and down. In summary, /5-decay rates of 25 Si, 33 Ar, 36-38 Ca, 
41 Ti, 59 Zn, 62 " 64 Ge, and 68 Se, become critical in Model K04-B1, whereas those of 25 Si, 68 Se, 
72 Kr, 76 Sr, 80 Zr, 92 < 94 Pd, and 101 - 104 Sn play a major role in Model K04-B2. As before, no 
noticeable effect remains when realistic uncertainties are adopted for these rates. 



3.1.5. Effect of the initial metallicity: Models K04-B3, KO4-B4, and K04-B5 

In the low metallicity Model K04-B3, the main nuclear path (X > 10~ 2 ) stops at 
96 Ru. In contrast, the higher metallicity models K04-B4 and K04-B5 reach 68 Ge and 72 Se, 
respectively. The most abundant species at the end of the bursts are: H (0.194), 4 He (0.021), 
68 Ge (0.205), 72 Se (0.132), 64 Zn (0.075), 76 Kr (0.073), 80 Sr (0.041), and 82 Sr (0.023), for Model 
K04-B3; 4 He (0.018), 60 Ni (0.696), 64 Zn (0.161), 56 Ni (0.051), 32 S (0.014), and 68 Ge (0.016), 
for Model K04-B4; and finally, 4 He (0.018), 56 Ni (0.273), 60 Ni (0.256), 39 K (0.057), 64 Zn 
(0.062), and 68 Ge (0.041), for Model K04-B5. Note that H has been seriously depleted in 
both high metallicity models, achieving mass fractions of 4 x 10 -14 and 10~ 15 , respectively. 

Concerning nuclear uncertainties, Model K04-B3 is characterized by 56 critical reactions, 
while Models K04-B4 and K04-B5, have 43 and 45, respectively. The most important ones 
for Model K04-B3 (Table 6) are 65 As(p,7) and 96 Ag(p,7), with a marginal role played by 
102 ' 103 In(p,7). For Model K04-B4 (Table 7), the most influential reactions are 12 C(a,7), 
30 S(o;,p), and to some extent, 30 P(a,p), and 65 As(p,7). Finally, for Model K04-B5 (Table 
8), the key reactions are 30 S(a,p), 25 Si(a,p), 59 Cu(p,7), 56 Ni(a,p), 29 S(a,p), 65 As(p,7), and 
12 C(o;,7), and to some extent, 61 Ga(p,7). Note that 15 0(a,7) is only important for Model 
K04-B3. The number of influential a-induced reactions seems to increase with metallicity: 
only 2 reactions of this type are important in Models K04 and K04-B3, whereas 22 and 17 
are important for Models K04-B4 and K04-B5, respectively. 

Finally, the most influential /3-decays from Model K04-B3 are 68 Se, 72 Kr, 76 Sr, 80 Zr, 
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88 Ru, 92 Pd, and "In. Model K04-B4 is characterized by the importance of the beta-decays 
of 33 > 34 Ar, 37 " 39 Ca, 42 Ti, 46 Cr, 49 Fe, 55 Ni, 58 Zn, 68 Se, 72 Kr, 76 Sr, and 80 Zr, whereas for Model 
K04-B5, the beta-decay rates of 21 Mg, 24 > 25 Si, 28 - 30 S, 33 > 34 Ar, 37 > 38 Ca, 64 Ge, and 68 Se, are 
the most important (see Table 11). Again, no effect shows up when realistic (ground-state) 
uncertainties are adopted for these rates. 



3.1.6. Effect of the peak temperature: Model K04-B6 vs. K04-B7 

In Model K04, the main nuclear path (X > 1CT 2 ) reached 96 Ru. Because of the lower 
peak temperature achieved in Model K04-B6, the main nuclear path reaches 82 Sr. In contrast, 
the higher temperatures achieved in Model K04-B7 drives the main path up to 103 Ag. 

The most abundant species at the end of the burst are now H (0.151), 4 He (0.034), 
64 Zn (0.375), 68 Ge (0.193), 60 Ni (0.051), 72 Se (0.074), 76 Kr (0.031), 80 Sr (0.015), and 82 Sr 
(0.011), for the lower peak temperature Model K04-B6 (Table 9); and H (0.460), 4 He (0.013), 
68 Ge (0.058), 72 Se (0.069), 76 Kr (0.048), 80 Sr (0.031), and 96 Ru (0.026), for Model K04-B7 
(Table 10). It is worth noting that Model K04-B7 ends with a larger amount of hydrogen 
than K04-B6 as a result of the major role played by photodisintegration reactions when the 
temperature exceeds ~ 2 GK. This is illustrated in Fig. 2, where the time evolution of the 
hydrogen mass-fraction is plotted for the two models K04-B6 (low T) and K04-B7 (high T). 
Indeed, the bump exhibited by Model K04-B7 during the 20 s after T peak (that is, when the 
temperature ranges between 1.5 and 2.5 GK) is caused by the protons released through (7, 
p) reactions on a number of nuclear species. Hence, the final hydrogen mass fraction may 
provide a diagnostic of the peak temperature achieved during the explosion (for bursts of 
similar duration). 

Concerning nuclear uncertainties, Model K04-B6 is characterized by 49 critical reactions, 
the most important ones being 61 Ga(p,7), 82 Zr(p,7), 65 As(p,7), 86 ' 87 Mo(p,7), 84 Nb(p,7), with 
a minor role played by 92 Ru(p,7) (Table 9). For Model K04-B7, we find 53 critical reactions, 
with the largest roles played by 69 Br(p,7), followed by 96 Ag(p,7), and 103 In(p,7) (Table 10). 
It is also worth mentioning that in both models only 2a — induced reactions appear to be 
influential (one being 15 0(a,7)). 

As expected, the most important /3-decay rates involve heavier species in Model K04-B7 
than in K04-B6, as a result of the former's larger peak temperature (see Table 11). Hence, 
the f3 -decay rates of 18 Ne, 64 Ge, 68 Se, 72 Kr, 76 Sr, and 80 Zr, are influential in Model K04-B6, 
whereas those of 68 Se, 72 Kr, 76 Sr, 80 Zr, 88 Ru, 92 Pd, 101 > 102 Sn, are important in Model K04-B7 
(again only for uncertainty factors of 10, up and down). 
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4. 



Monte Carlo Simulations 



As discussed earlier, realistic simulations of Type I X-ray bursts have revealed a dramatic 
extension of the nuclear path, reaching the proton-drip line at many points, and extending 
typically up to A=50-60 (Fisker et al. 2008), or eventually, up to the SnSbTe-mass region 
(Schatz et al. 2001), or beyond (Koike et al. 2004). The overall number of nuclear processes 
is huge, since many different reactions are involved in the synthesis and/or destruction of 
any particular isotope during the explosion. 

Recently, it has been claimed (Roberts et al. 2006) that, because of the coupling of so 
many different channels in XRBs, traditional sensitivity studies (like the one presented in 
Section 3), in which only one reaction is varied while the others are held constant, cannot 
properly address all the important correlations between rate uncertainties and nucleosyn- 
thetic predictions, leading to wrong (or at least, biased) conclusions. 

Here, we will examine the impact of simultaneously varying all reaction rates in the 
highly-coupled environment characteristic of an XRB through Monte Carlo simulations, and 
will discuss the feasibility of traditional sensitivity studies through a direct comparison. 



The simultaneous variation of reaction rates relies on the generation of a set of pseudo- 
random numbers that follow a given probability density function. 

Since nuclear reaction rates are positive quantities, we assume that the random enhance- 
ment factors that are applied to each individual reaction follow a log-normal probability 
function (note that other distributions, such as Gaussian, may give negative values). The 
log-normal probability function adopted has the form: 



where fi and a are the mean and the standard deviation of ln(x), respectively^. 

This log-normal distribution is implemented in such a way that the geometric mean of 
all enhancement factors is fi* — 1, and the probability to generate a random number between 



3 By definition, if x is log-normally distributed, ln(x) is normally distributed. Note that /i* = e' 1 is the 
geometric mean (as well as the median) of the log-normal distribution, whereas a* = e a is the geometric 
standard deviation. 



4.1. 



The Monte Carlo code 
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0.1 and 10 (that is, in the range [fi* / (a*) 2 , fj,* x (a*) 2 ]), to match the factors adopted for 
our sensitivity study in Section 3, is chosen to be 95.5%. This is accomplished by setting 
a* = V4~0. From the many algorithms available in the literature for generating pseudo- 
random numbers, we have used that from Guardiola et al. (1995). 

The set of enhancement factors is applied to all nuclear reactions included in the network; 
the same enhancement factor was applied to a forward-reverse reaction paii0. This is coupled 
to a post-processing nucleosynthesis code that allows us to compute the final yields for a 
particular temperature and density versus time profile, in a similar way as described in 
Sections 2 & 3. The process is then iterated for a pre-selected number of trials, in order to 
achieve statistically sound results. The series of computations performed in this work rely on 
1000 trials. Although previous Monte Carlo simulations have been based on a larger number 
of trials (10,000, in Smith et al. 2002 and Hix et al. 2003; 50,000, in Roberts et al. 2006), 
we feel our choice is sufficient for the goals of this manuscript (see sections 4.2 and 5). 



4.2. Results 

We have applied the Monte Carlo technique to our set of models discussed in Section 
3. However, since results are qualitatively similar in all cases we will restrict our discussion 
here to Model K04 (Koike et al. 2004). 

Before a thorough analysis of the results is made, several aspects must be considered 
in detail. First, the use of realistic uncertainties affecting all the rates is critical in Monte 
Carlo simulations, since the simultaneous variation of all reactions may otherwise result in 
inflated uncertainties in the final yields (Smith et al. 2002). Because of this, studies based on 
individual reaction-rate variations are more easily interpreted. This is illustrated in Figures 
3 & 4, where relative abundances and their normalized geometric standard deviations (see 
Limpert et al. 2001, and Carobbi et al. 2003) are plotted as a function of the mass number, 
for Model K04. The relative abundances correspond to the geometric mean abundances 
obtained from the set of Monte Carlo trials normalized to abundances found with standard 
rates. In turn, the normalized geometric standard deviations define the interval for which the 
probability to enclose any possible geometric mean abundance for a given species matches a 
certain value (95.5%, in this study). 



A rigorous Monte Carlo procedure would require a random sampling of the reaction Q-value for the 
calculation of the reverse reaction rate. We did not address this extra complication, which we leave to future 
investigations. 
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The simultaneous variation of all nuclear processes in bulk between 0.1 and 10, according 
to our log-normal distribution, results in large uncertainty bars for many isotopes (Figure 
3). Notice, however, that these final abundance uncertainties are overestimated since all 
/3 + -decays, as well as some important, relatively well-known reactions, such as the triple-a, 
were allowed to vary in this way. As discussed in Section 3, more realistic uncertainties 
must be used instead, since these rates are known with better precision (usually better than 
~ 30%, for /3-decay rates). Indeed, when all nuclear processes, except the 3a and all the 
/9-decay rates, are allowed to vary between roughly 0.1 and 10 in the Monte Carlo study, 
the overall uncertainties in the final yields decrease dramatically (Figure 4). Similar results 
have been achieved for Model F08 (see Figures 5 & 6). 

A second warning associated with Monte Carlo simulations involves the reduced subset 
of reactions whose variation affects the overall energy production (see Section 5). In individ- 
ual reaction-rate variation studies, these reactions can be appropiately flagged for separate, 
detailed analysis with better numerical tools (semi-analytical or hydrodynamic codes that 
can properly address changes in the temperature and density profiles driven by variations 
in the total energy output). In Monte Carlo studies, however, one cannot simply remove 
those trials in which the overall energy production is modified as this would affect the input 
distribution of enhancement factors (which are assumed to be random). Therefore, while re- 
sults from individual reaction-rate variations are not corrupted by these effects, Monte Carlo 
simulations cannot disentangle this from the overall analysis and hence the interpretation of 
the results has to be taken with caution. 

Figure 4 shows the impact of the simultaneous variation of all rates (except for the 
3a and all the /3-decay rates) on the final yields, for the temperature and density versus 
time profiles of Model K04 (Koike et al. 2004). It illustrates the interplay of multiple 
nuclear processes in the highly coupled environment of an XRB. The identification of the key 
reactions, whose uncertainties have the largest impact on the final yields, is more complicated 
than in traditional sensitivity studies. Figure 4 shows indeed which isotopes are mostly 
affected by the uncertainties associated with the rates. However, the identification of those 
specific reactions that are perhaps most responsible for those changes is, by no means, 
straightforward. Following Smith et al. (2002), Hix et al. (2003), and Roberts et al. (2006), 
we have searched for possible correlations between variations in the final abundance of a 
specific nucleus and each nuclear reaction rate that was varied in the Monte Carlo routine. 

In many cases, correlations are small, as illustrated in Figure 7: although the relative 
abundance of 69 Ge shows a significant variability, it cannot be attributed to the uncertainty 
affecting the 56 Ni(a, p) 59 Cu (and 59 Cu(p,a) 56 Ni) rate, since the fit has a correlation coeffi- 
cient of only r=0.101. In sharp contrast, Figs. 8 & 9 show a large correlation between the 
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69 Ge abundance and 69 Se(p, 7) 70 Br, as well as between 72 Se and 65 As(p, 7) 66 Se. However, 
despite the large correlation, varying the 65 As(p, 7) rate by a factor of 10 has no dramatic 
impact on the final 72 Se yield (Fig. 9). Hence, a key reaction in the Monte Carlo approach 
must have a significant correlation with the yield of an isotope (as indicated by the correla- 
tion coefficient) as well as a pronounced impact on the final abundance of an isotope as its 
rate is varied. 

The most important correlations from the Monte Carlo study applied to Model K04, as 
characterized by fits with correlation coefficients r > 0.5, are listed in Table 15. From the 
50 reactions listed in the table, all but 5 (namely, 34 Ar(p, 7), 53 Co(p, 7), 62 Ga(p, 7), 69 Br(p, 
7), and 80 Y(p, 7)) were previously identified in the study based on individual reaction-rate 
variations (Table 1). 

Table 16 displays reactions that have the largest effects on the final yields for Model 
K04, according to the individual reaction-rate variation study (column 2) and Monte Carlo 
simulations (column 3). Reactions included in Table 16 were restricted to those that, in either 
the individual reaction-rate variation study or in the Monte Carlo simulations, affected the 
final abundances of the isotopes by at least a factor of 2 (for the former, see Table 1; for 
the latter, reactions from Table 15 without brackets in the slope column). Qualitatively, 
results are in excellent agreement, except at the very end of the network (nuclei with A > 
97). Indeed, most of the discrepancies correspond to reactions that were indeed identified in 
Table 15, but that affect the final abundances by less than a factor of 2 when enhancement 
factors ranging between 0.1 and 10 are considered. These discrepancies might be related 
to the correlations that the individual reaction-rate variation study are not capable to deal 
with. 

Similar agreement is reported on the comparison of both techniques when applied to 
Model F08 (Fisker et al. 2008), as illustrated in Tables 17 & 18. 

Finally, we would like to stress that our preliminary Monte Carlo studies rely on 1,000 
trials, a limited number that calls for a justification. To this end, we have extended the 
study of Model K04 to 4,200 trials. Indeed, no significant differences, neither in the list of 
key reactions, nor on the correlation coefficients and slopes have been found. 

5. Discussion and Conclusions 

Table 19 summarizes the most important reactions collected from Tables 1-10. For 
the sake of brevity we have restricted Table 19 and the dicussion here to those reactions 
that affect the yield of 3 species or more in any of our 10 models, in the individual- variation 
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studies. We have also carefully identified reactions that were seen to modify the XRB energy 
output when their rates were varied by a factor of 10, up or down. This is explicitly indicated 
in Table 20, which lists the subset of reactions with any impact on XRB yields (reactions 
identified in Tables 1 to 10, that affect, at least, one isotope) that simultaneously modify the 
overall energy output by more than 5% at some point during the burst, when their nominal 
rates are varied by a factor of 10, up or down. This table has to be taken as a warning of 
the limitations of post-processing techniques (both for individual- variation and Monte Carlo 
studies). Furthermore, we have identified some additional reactions in our studies which 
affect the energy output, but remarkably did not affect any yields in any of our models (for 
instance, 14 0(a, p), 27 Si(p, 7), 31 S(p, 7), or 35 S(p, 7)). Several aspects are worth noting 
here. First, the total number of reactions affecting the energy output, for any model, is 
small. Second, as indicated in Table 20, some of those reactions are known with better 
precision than a factor of 10. Third, there is no way to overcome this problem in the context 
of post-processing calculations. Indeed, a self-consistent analysis with a hydrodynamic code 
capable of self-adjusting both the temperature and the density of the stellar envelope seems 
mandatory to address this issue, for the few reactions of concern. And fourth, the presence of 
reactions affecting the energy output is particularly dangerous within a Monte Carlo context, 
as results rely on the simultaneous variation of all rates, whereas most of the results obtained 
from an individual reaction-rate variation approach remain unaffected. Note, however that 
Tables 19 & 20 reveal a very limited number of reactions that were seen to affect the energy 
output for Models K04 and F08, giving support to the results obtained in our Monte-Carlo 
studies. 

It is also worth mentioning that no discussion involving the triple-a reaction or any beta- 
decay rate has been made here, since these reactions, when varied within realistic uncertainty 
limits, have no effect, neither on yields, nor on the nuclear energy output. The same applies 
to Q-value variations, hence results listed in Tables 13 & 14 are not affected by variations of 
the energy output. 

5.1. Overview of the most influential reactions 

The main purpose of this Section is to elaborate on the reactions listed in Tables 19 & 
20, namely, the most influential reactions found in this study. 

First, we will focus on those rates that were drawn from theoretical estimates, due to 
insufficient or unavailable experimental information. Their main characteristics are summa- 
rized in Table 21. To help illustrate systematic uncertainities associated with these important 
reactions, we have compared the theoretical rates adopted in our network with those found 



-22 - 



in the recent REACLIBvO compilation.] over the temperature range covered in our studies. 
Overall, the agreement is quite good; however, we note that even when using the same ba- 
sic Hauser-Feshbach code (NON-SMOKER) along with similar Q-values, a difference in the 
rates as large as a factor of ~3 is obtained. Differences as large as a factor of ~4 are seen 
between rates from different statistical model codes, but using similar Q-values. Moreover, 
the factor of 10 disagreement for the 61 Ga(p, 7) rate arises from the comparison between a 
NON-SMOKER result and a shell-model calculation. The magnitude of these discrepancies 
lends support to our choice of varying rates by a factor of 10 rather than by a significantly 
larger factor. 

We will now focus on the rates listed in Tables 19 & 20 that have been determined 
experimentally, and we will assess whether a factor of 10 variation is reasonable for these 
rates. For those cases where a smaller uncertainty is justified, we have performed addi- 
tional post-processing calculations to supplement results from Tables 1-10 and 20 (namely, 
to determine the impact of individual rate variations by smaller factors on yields and on the 
overall nuclear energy outuput). For reference, we will continue to compare the rate adopted 
in this work with that used in the REACLIBvO compilation. Although we do not discuss ex- 
perimental information (if available) for other reactions listed in Tables 1-10 (namely, those 
that affected less than 3 isotopes in any model), this must of course be considered by anyone 
examining reactions beyond the most influential ones listed in Table 19. 

• 12 C(a, 7) 16 0: We have used the Kunz et al. (2002) experimental rate which, over 
the range of temperatures used in our models, agrees to better than 20% with the 
recommended NACRE (Angulo et al. 1999) rate, and to better than a factor of ~2 
with the recommended Buchmann (1996) rate adopted in REACLIBvO. A factor of 
~3 variation in our rate would cover the limits given in both NACRE and Buchmann 
(1996). Hence, we have varied this rate individually by a factor of 3, up and down, 
and tested its effect for those models in which a factor of 10 variation, as discussed 
in Section 3, had an impact (Models F08, K04-B2, K04-B4, and K04-B5). Indeed, 
only Model K04-B2 (where the 12 C yield is affected by a factor 0.47 when this rate is 
multiplied by 3) and K04-B5 (where 20 Ne and 24 Mg are affected by factors 0.44 and 
0.47, respectively, when the rate is reduced by a factor 3) reveal changes in the final 
yields (note that variation of this reaction by a factor of 10 did not affect the overall 
nuclear energy output in any of our models). 

• 18 Ne(a, p) 21 Na: We use the Chen et al. (2001) experimental rate, which agrees to 
a factor ~3 with the Gorres et al. (1995) Hauser-Feshbach (SMOKER) calculation 



5 http:/ /www. nscl.msu.edu/ nero/db/ 
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adopted in REACLIBvO. There is, however, additional data from Groombridge et al. 
(2002). Indeed, using the information from Chen et al. (2001), for E r < 1.7 MeV, and 
that from Groombridge et al. (2002), for E r > 1.98 MeV, we find a rate which deviates 
from the Chen et al. rate by 30%, a factor of ~3 and a factor of ~7 at 1.0, 1.4, and 2 
GK, respectively. Since we found this reaction to affect yields and/or the total nuclear 
energy in Models K04-B1 and K04-B6, neither of which reach temperatures above 1.36 
GK, we have restricted our analysis to a variation of the Chen et al. rate by a factor 
of 3, up and down. No impact on the yields is found when varying the rate as such; 
however, Model K04-B1 shows some variation of the total energy output when the rate 
is multiplied by 3. 

• 15 0(a, 7) 19 Ne: We use the Davids et al. (2003) rate, which agrees to factor of ~3 
with the Hahn et al. (1996) rate used in REACLIBvO, over the temperatures spanned 
by our models. If we combine the new information reported by Tan et al. (2007), for 
E x < 4.55 MeV, with the information in Davids et al. (2003), for the states at E x = 
4.600, 4.712 and 5.092 MeV, we calculate a new rate that agrees to a factor ~2 with 
the original Davids et al. rate, over the relevant temperatures. Varying our rate within 
a factor ~3 would cover the uncertainty limits in the Tan-Davids rate calculation. We 
have tested the impact of varying our l5 0(a, 7) rate by a factor 3 up and down, for 
Models K04, K04-B1, and K04-B6. We find that the total nuclear energy is affected 
only in the early stages of Model K04, when the rate is multiplied by 3. Concerning 
the yields, only 15 N is affected (consistent with Tables 1, 4 and 9): when the rate is 
reduced by a factor 3 the 15 N normalized yields in Models K04, K04-B1 and K04-B6 
are 3.13, 2.75, and 3.17, respectively; when the rate is increased by a factor 3, no effect 
is seen except in Model K04-B1, where the normalized yield of 15 N decreases to 0.35. 

• 23 Al(p, 7 ) 24 Si: We use the Schatz et al. (1997) rate, as does REACLIBvO. This rate 
is based on both theoretical estimates and measurements of excited states in 24 Si. Its 
uncertainty spans up to 3 orders of magnitude, for typical XRB temperatures. For this 
reason, we deem our results from varying this rate by a factor 10 to be adequate. 

• 24 Mg(a, p) 27 Al: Both REACLIBvO and our network rely on the experimental rate 
reported in Iliadis et al. (2001). Above 0.3 GK, this rate is reported to be uncertain 
by only ±20%. No effect on the yields is found when the rate is varied by 20% in our 
post-processing calculations (Model K04-B2; see Table 20). The total nuclear energy 
is, however, affected when this rate is increased by 20%, in that particular model. 

• 26s Al(p, 7) 27 Si: We use the (unpublished) rate from the PhD thesis of Vogelaar (1989), 
as does REACLIBvO (the new TRIUMF measurement -Ruiz et al. (2006)- on the 
184 keV resonance does not affect the rate at typical XRB temperatures). According 
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to NACRE, the experimental information is only sufficient to determine the rate for 
T < 0.9 GK (they use Hauser-Feshbach calculations to extend this rate to higher 
temperatures). Our rate agrees with the NACRE recommended rate to 30% over the 
temperature range covered by our models; varying our rate by a factor of 2 up and 
down would cover the uncertainties reported in NACRE. Assuming such a degree of 
uncertainty in our rate for Model F08 (see Table 20) leads to no effect on any yield, 
nor on the total nuclear energy. 

• 28 Si(o;, p) 31 P: The rate used in our network and that of REACLIBvO is from the 
compilation of Iliadis et al. (2001). An uncertainty of ±20% is assigned to this rate 
above 0.2 GK. Actually, varying this rate by ± 20% in Model K04-B4 (see Table 20) 
has no effect on any yields, nor on the total nuclear energy. 

• 32 S(a, p) 35 Cl: We use the Iliadis et al. (2001) experimental rate, as does REACLIBvO. 
The uncertainty of this rate, over XRB typical temperatures, spans up to 3 orders of 
magnitude (Iliadis et al. 2001). Consequently, we feel that our results, based on a 
variation of this rate by a factor of 10 (see Table 20), are reliable. 

• 35 Cl(p, 7) 36 Ar: Again, we adopt the Iliadis et al. (2001) experimental rate, as does 
REACLIBvO. The uncertainty in this rate is ~20% above 0.2 GK (Iliadis et al. 2001). 
We find no effect on any yield when varying this rate by ±20% in our post-processing 
calculations for Model K04-B2 (see Table 20). However, we remarkably find that 
increasing this rate by 20% does indeed affect the total nuclear energy at the early 
stages of the TNR in this model. 

Finally, because of past interest in the literature, we discuss some specific reactions that 
we found to affect the total nuclear energy when their rates were varied individually by a 
factor of 10, but affected no yields whatsoever (see Section 5.1): 31 S(p, 7), 35 Ar(p, 7), 27 Si(p, 
7), and 14 0(a, p). 

• 31 S(p, 7) 32 C1, 35 Ar(p, 7) 36 K, and 27 Si(p, 7) 28 P: For these 3 reactions, we have used 
the rates reported in Iliadis et al. (1999), determined through the use of experimental 
information when available (e.g., excitation energies), along with calculations and in- 
formation from the respective mirror nuclei. Judging from the uncertainties presented 
in that work, we conclude that uncertainty factors of 2, 3 and 2 (up and down) are more 
reasonable for the 31 S(p, 7), 35 Ar(p, 7), and 27 Si(p, 7) rates over XRB temperatures, 
respectively. Accordingly, we individually varied each of these reactions within those 
limits for the models affected (i.e., Models K04-B1 and K04-B7 for 31 S(p, 7); Models 
K04-B2 and K04-B7 for 35 Ar(p, 7); and Models K04-B4 and K04-B7 for 27 Si(p, 7)) to 
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determine the resulting impact on the total nuclear energy, as no yields were affected 
by varying any of these by a factor of 10. Increasing the 31 S(p, 7) rate by a factor 
of 2 changed the total nuclear energy (by at least 5% at some point of the burst) in 
Model K04-B1; increasing the 35 Ar(p, 7) rate by a factor 3 changed the nuclear energy 
in Model K04-B2; and increasing the 27 Si(p, 7) rate by a factor 2 changed the nuclear 
energy in Model K04-B7. The impact of these reactions on XRB light curves has been 
examined in Iliadis et al. (1999) and in Thielemann et al. (2001). 

• 14 0(a, p) 17 F: This reaction is of critical importance for breakout from the hot CNO- 
cycle. We use the rate from Blackmon et al. (2003), which is larger than the Hahn et al. 
(1996) rate, adopted by REACLIBvO, by a factor of ~10 at typical XRB temperatures. 
This increased rate is due to the inclusion of some of the contributions from the 17 F* 
exit channel by Blackmon et al.; branches for states of E x > 7 MeV in 18 Ne were 
not measured though. Actually, Notani et al. (2004a, b) observed what could be the 
population of 17 F* through a state at E x ( 18 Ne) ~ 7.1 MeV, but the interpretation of 
their results has been questioned by Fu et al. (2007). As the studies of Blackmon et 
al. and Notani et al. are both published in only preliminary forms, and taking into 
account the argument of Fu et al., we find it difficult to evaluate the uncertainty in 
this rate. Variations by a factor of 10 affected nuclear energy in 5 of our 10 models: 
K04, F08, K04-B1, K04-B3 and K04-B4. Further efforts to constrain this rate, based 
on analyses of previous measurements and/or new measurements, would certainly be 
desirable. 

In summary, we have identified a very limited number of reactions (see Tables 19-21 
and this section) that play a significant role in our XRB nucleosynthesis studies. Indeed, 
our results can help to guide and motivate future measurements by experimental nuclear 
physicists. Stellar modelers, as well, may tackle the challenge to properly address the role 
played by the few reactions flagged as affecting the overall energy output, an aspect that lies 
beyond the possibilities offered by post-processing calculations. 

Finally, we have compared the two general approaches to the sensitivity study: indi- 
vidual reaction-rate variations, and simultaneous variation of all rates through Monte Carlo 
techniques. The highly coupled environment characteristic of an XRB provides the oppor- 
tunity to test for the effects of correlations in the uncertainties of reaction rates. We find 
similar results from both approaches, with minor differences attributed to such correlation 
effects. 

This work has been partially supported by the Spanish MEC grant AYA2007-66256, 
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Table 1. Final abundance ratios, X { j X ijStd , for Model K04 (Koike et al. 2004), resulting 
from reaction-rate variations by a factor of 10 (up and down). 



Reaction 


Isotope 


10 


0.1 


15 0(a, 7 ) 19 Ne 


15 N 




8.30 


18 Ne(a, p) 21 Na 


15 N 


0.43 






18p 


0.43 




31 Cl(p, 7 ) 32 Ar 


30 Si 


0.44 




43 V(p, 7 ) 44 Cr 


44 Ti 




0.49 


47 Mn(p, 7 ) 48 Fe 


46 Ti 


0.31 




56 Ni(p, 7 ) 57 Cu 


56 Ni 




7.78 


56 Cu(p, 7 ) 57 Zn 


55 Co 




2.64 


57 Cu(p, 7 ) 58 Zn 


57 Ni 




2.79 


58 Cu(p, 7 ) 59 Zn 


58 Ni 




2.05 


59 Cu(p, 7 ) 60 Zn 


59 Ni 




2.02 


61 Zn(p, 7 ) 62 Ga 


61 Cu 




2.35 


61 Ga(p, 7 ) 62 Ge 


eo Ni 


0.16 


4.10 




61 Cu 




4.25 


63 Ga(p, 7 ) 64 Ge 


63 Cu 




3.29 


65 Ge(p, 7 ) 66 As 


65 Zn 




8.95 


65 As(p, 7 ) 66 Se 


64 Zn 


0.35 


3.50 




65 Zn 


0.37 


2.93 




66 Ge 


0.44 


2.77 




67 Ga 


0.40 


3.22 




81 Rb 




0.50 




88 Zr 




0.49 




89 Nb 




0.49 




90 Mo 




0.48 




91 Nb 




0.47 




92 Mo 




0.47 




93 Tc 




0.46 




94 Mo 




0.45 




95 Ru 




0.44 




96 Ru 




0.43 




97 Ru 




0.40 




98 Ru 




0.39 




"Rh 




0.38 
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Table 1 — Continued 



Reaction Isotope 10 0.1 





lOOp^ 




0.38 




ioi pd 




0.38 




102p d 




0.37 




103 Ag 




0.38 




104 Ag 




0.40 


66 Ge(p, 7 ) 67 As 


66 Ge 




4.86 


66 As(p, 7 ) 67 Se 


66 Ge 


0.18 




67 As(p, 7 ) 68 Se 


67 Ga 


0.47 


5.71 


69 Se(p, 7 ) 70 Br 


69 Ge 


0.10 


8.21 


70 Se(p, 7 ) 71 Br 


70 Ge 




5.18 


70 Br(p, 7 ) 71 Kr 


70 Ge 


0.28 




71 Br(p, 7 ) 72 Kr 


71 As 


0.20 


7.11 


73 Kr(p, 7 ) 74 Rb 


73 Se 


0.10 


7.50 


74 Kr(p, 7 ) 75 Rb 


74 Se 


0.30 


6.51 


74 Rb(p, 7 ) 75 Sr 


74 Se 


0.40 




75 Rb(p, 7 ) 76 Sr 


75 Br 


0.15 


6.35 


77 Sr(p, 7 ) 78 Y 


77 Kr 


0.11 


5.72 


78 Sr(p, 7 ) 79 Y 


78 Kr 


0.17 


5.82 


79 Y(p, 7 ) 80 Zr 


79 Kr 


0.12 


5.45 


81 Zr(p, 7 ) 82 Nb 


81 Rb 


0.15 


3.15 


82 Zr(p, 7 ) 83 Nb 


82 Sr 


0.17 


2.92 


83 Nb(p, 7 ) 84 Mo 


83 Sr 


0.27 




84 Nb(p, 7 ) 85 Mo 


84 Sr 


0.20 


2.44 


85 Nb(p, 7 ) 86 Mo 


85y 


0.43 




85 Mo(p, 7 ) 86 Tc 


85y 


0.48 




86 Mo(p, 7 ) 87 Tc 


86 Zr 


0.23 


2.43 


87 Mo(p, 7 ) 88 Tc 


87 Zr 


0.32 




88 Tc(p, 7 ) 89 Ru 


88 Zr 


0.37 


2.04 


89 Tc(p, 7 ) 90 Ru 


89 Nb 


0.24 


2.15 


90 Ru(p, 7 ) 91 Rh 


90 Mo 


0.35 




91 Ru(p, 7 ) 92 Rh 


91 Nb 


0.37 




92 Ru(p, 7 ) 93 Rh 


92 Mo 


0.32 




93 Rh(p, 7 ) 94 Pd 


93 Tc 


0.24 


2.05 
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Table 1 — Continued 



Reaction 


Isotope 


10 


0.1 


94 Rh(p, 7 ) 95 Pd 


94 Mo 


0.47 




94 Pd(p, 7 ) 95 Ag 


94 Mo 


0.45 




95 Pd(p, 7 ) 96 Ag 


95 Ru 


0.31 




96 Ag(p, 7 ) 97 Cd 


96 Ru 


0.48 






97 Ru 




0.41 








U.oo 




"Rh 




0.43 




ioo pd 




0.48 


97 Ag(p, 7 ) 98 Cd 


97 Ru 


0.26 


2.13 


98 Ag(p, 7 )"Cd 


98 Ru 


0.45 




98 Cd(p, 7 )"ln 


98 Ru 


0.43 




"Cd(p, 7 ) 100 In 


99 Rh 


0.32 




100 Cd(p, 7 ) 101 In 


100 Pd 


0.43 




101 In(p, 7 ) 102 Sn 


ioi pd 


0.34 




102 In(p, 7 ) 103 Sn 


102p d 


0.36 






103 Ag 




0.40 




104 Ag 




0.48 


103 In(p, 7 ) 104 Sn 


104 Ag 


2.80 


0.24 




105 Ag 


2.00 




104 In(p, 7 ) 105 Sn 


105 Ag 


4.13 





-35 - 



Table 2. Same as Table 1, for Model SOI (Schatz el at. 2001). 



rvcclCLlUIl 


1 C f~\~i~ f~\T~\ O 


1 n 

1U 


n 1 

U. 1 


15 0(a, 7 ) 19 Ne 


15 N 




8.84 


52 Fe(p, 7 ) 53 Co 


52 Fe 




8.48 


52 Co(p, 7 ) 53 Ni 

\I 7 1 / 


52 Fe 




2.03 


56 Ni(p, 7 ) 57 Cu 


56 Ni 




8.33 


56 Ni(a, p) 59 Cu 




0.36 






52 Fe 


2.40 






56 Ni 


2.31 






69 Ge 


2.27 






70 Ge 


2.19 






71 As 


2.02 






73 Se 


2.25 






74 Se 


2.19 




57 Ni(p, 7 ) 58 Cu 


57 Ni 




2.50 


57 Cu(p, 7 ) 58 Zn 


57 Ni 




2.09 


58 Cu(p, 7 ) 59 Zn 


58 Ni 




4.92 


59 Cu(p, 7 ) 60 Zn 


X H 




0.35 




52 Fe 




2.45 




56 Ni 




2.35 




59 Ni 




6.60 




69 Ge 




2.32 




70 Ge 




2.24 




71 As 




2.07 




73 Se 




2.29 




74 Se 




2.24 


61 Zn(p, 7 ) 62 Ga 


61 Cu 




5.59 


61 Ga(p, 7 ) 62 Ge 


eo Ni 


0.29 




62 Zn(p, 7 ) 63 Ga 


62 Zn 




3.52 


63 Ga(p, 7 ) 64 Ge 


63 Cu 


0.31 


6.80 


65 Ge(p, 7 ) 66 As 


65 Zn 




7.90 


66 Ge(p, 7 ) 67 As 


66 Ge 


0.22 


6.19 


66 As(p, 7 ) 67 Se 


66 Ge 


0.40 




67 As(p, 7 ) 68 Se 


67 Ga 


0.13 


6.71 


69 Se(p, 7 ) 70 Br 


69 Ge 


0.11 


6.14 
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Table 2 — Continued 



Reaction Isotope 10 0.1 



69 Br(p, 7 ) 70 Kr 


72 Se 




2.06 




76 Kr 




2.06 


70 Se(p, 7 ) 71 Br 


70 Ge 


0.21 


7.01 


71 Br(p, 7 ) 72 Kr 


71 As 


0.11 


5.53 


73 Kr(p, 7 ) 74 Rb 


73 Se 


0.11 


5.53 


74 Kr(p, 7 ) 75 Rb 


74 Se 


0.14 


6.84 


75 Rb(p, 7 ) 76 Sr 


75 Br 


0.12 


4.54 


77 Sr(p, 7 ) 78 Y 


77 Kr 


0.16 


3.17 


78 Sr(p, 7 ) 79 Y 


78 Kr 


0.12 


5.42 


79 Y(p, 7 ) 80 Zr 


79 Kr 


0.14 


3.43 


80 Y(p, 7 ) 81 Zr 


80 Sr 


0.43 


2.84 


81 Y(p, 7 ) 82 Zr 


81 Rb 




2.15 


81 Zr(p, 7 ) 82 Nb 


81 Rb 


0.27 




82 Zr(p, 7 ) 83 Nb 


82 Sr 


0.20 


2.25 


83 Zr p, 7 84 Nb 


83 Sr 


0.37 




84 Zr(p, 7 ) 85 Nb 


84 Sr 


0.46 




84 Nb(p, 7 ) 85 Mo 


84 Sr 


0.42 




85 Nb(p, 7 ) 86 Mo 


85y 


0.29 




86 Mo(p, 7 ) 87 Tc 


86 Zr 


0.40 




87 Mo(p, 7 ) 88 Tc 


87 Zr 


0.36 




88 Mo(p, 7 ) 89 Tc 


88 Zr 


0.27 




89 Tc(p, 7 ) 90 Ru 


89 Nb 


0.32 




90 Mo(p, 7 ) 91 Tc 


90 Mo 


0.45 




91 Tc(p, 7 ) 92 Ru 


91 Nb 


0.30 




92 Ru(p, 7 ) 93 Rh 


92 Mo 


0.23 




93 Ru(p, 7 ) 94 Rh 


93 Tc 


0.47 




93 Rh(p, 7 ) 94 Pd 


93 Tc 


0.30 


2.29 


94 Rh(p, 7 ) 95 Pd 


94 Mo 


0.28 




94 Pd(p, 7 ) 95 Ag 


94 Mo 


0.49 




95 Rh(p, 7 ) 96 Pd 


95 Ru 


0.46 




95 Pd(p, 7 ) 96 Ag 


95 Ru 


0.32 




96 Ag(p, 7 ) 97 Cd 


96 Ru 


0.46 


2.29 


97 Ag(p, 7 ) 98 Cd 


97 Ru 


0.22 





-37- 



Table 2 — Continued 



Reaction Isotope 10 0.1 



97 Cd(p, 7 ) 98 In 


97 Ru 


0.43 




98 Ag(p, 7 )"Cd 


98 Ru 


0.42 




98 Cd(p, 7 )"ln 


98 Ru 


0.40 


2 


"Cd(p, 7 ) 100 In 


"Rh 


0.28 




100 Cd(p, 7 ) 101 In 


ioo pd 


0.49 




100 In(p, 7 ) 101 Sn 


ioo pd 


0.44 




101 In(p, 7 ) 102 Sn 


ioi pd 


0.30 




102 In(p, 7 ) 103 Sn 


102p d 


0.36 




103 In(p, 7 ) 104 Sn 


103 Ag 


0.46 




103 Sn(a, p) 106 Sb 


106 Cd 


0.34 






107 Cd 


0.24 






108 Cd 


0.21 






109 In 


0.21 




105 Sn(p, 7 ) 106 Sb 


105 Ag 


0.33 




107 Sn(p, 7 ) 108 Sb 


108 Cd 


2.32 




108 Sb(p, 7 ) 109 Te 


109j n 


3.33 




109 Sb(p, 7 ) 110 Te 


109 In 


0.46 


2 
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Table 3. Same as Table 1, for Model F08 (Fisker el at. 2008). 



Reaction Isotope 10 0.1 



12 C(p, 7) 13 N 


12 C 


0.41 




12 C(a -y) 16 


12 C 


0.42 






16Q 


3.27 


0.36 




20 Ne 


3.03 


0.42 




24 M K 


2.40 




m O(a V) 20 Ne 

V J I KJL , J J INC 


16Q 


1 1 


8 36 




20 Ne 


0.10 


9.84 


22 Me(a d) 25 A1 


25 Me 




8.81 




27 A1 




2.23 


24 Mefo <-y) 25 Al 


25 Me 


3.02 


0.43 


24 Me(a 7) 28 Si 


24 Me 


0.18 


4.16 




25 Me 


0.24 


4.10 


25 A1(« n^Si 

ill! Lt^ ) kJ 1 


25 Me 


0.22 


2 36 


26fl AlfD V) 27 Si 


27 A1 


2.71 


0.40 


26 sAl(a, p) 29 Si 


26 Mg 




2.34 




26 9A1 


0.08 


2.38 




27 A1 


0.40 




26 Si(a p) 29 P 


26g Al 


0.08 






27 A1 


0.23 




27 Al(a, p) 30 Si 


27 A1 


0.32 




27 Si(ct, p) 30 P 


27 A1 


0.30 




28 Si(a, 7) 32 S 


28 Si 


0.43 






32 S 


3.73 


0.37 


29 Si(a, 7 ) 33 S 


33 S 


2.15 




30 P(p, 7 ) 31 S 


31p 


2.38 


0.29 


30 P(«, p) 33 S 


33 s 


2.29 




34 Cl(p, 7 ) 35 Ar 


35Q 


2.07 




36 Ar(p, 7 ) 37 K 


37 Ar 


0.40 


2.76 


39 K(p, 7 ) 40 Ca 


40 Ca 


2.15 




44 Ti(p, 7 ) 45 V 


45 Ti 


0.18 


2.32 


44 V(p, 7 ) 45 Cr 


45 Ti 




0.47 


45 V(p, 7 ) 46 Cr 


45 Ti 




2.25 


46 V(p, 7 ) 47 Cr 


47 Ti 


2.22 
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Table 3 — Continued 



Reaction Isotope 10 0.1 



47 Cr(p, 7 ) 48 Mn 


47 Ti 


0.42 


2.78 


49 Mn(p, 7 ) 50 Fe 


49y 




2.19 


52 Fe(p, 7 ) 53 Co 


52p e 


0.31 


3.06 




53 Mn 




0.48 


53 Co(p, 7 ) 54 Ni 


53 Mn 


0.18 


4.06 


54 Co(p, 7) 55 Ni 


54p e 




2.01 


57 Ni(p, 7 ) 58 Cu 


57 Ni 




7.71 


57 Cu(p, 7 ) 58 Zn 


56 Ni 




2.07 




57 Ni 


0.46 


2.31 




58 Ni 


3.31 


0.37 




59 Ni 


2.20 


0.47 


58 Cu(p, 7 ) 59 Zn 


58 Ni 


0.44 


5.57 


59 Cu(p, 7 ) 60 Zn 


59 Ni 




2.88 


61 Zn(p, 7 ) 62 Ga 


61 Cu 


0.21 


3.59 




62 Zn 




0.49 


61 Ga(p, 7 ) 62 Ge 


64 Zn 


3.11 






65 Zn 


3.06 






66 Ge 


3.06 






67 Ga 


3.13 




62 Ga(p, 7 ) 63 Ge 


62 Zn 


0.23 


2.39 


63 Ga(p, 7 ) 64 Ge 


63 Cu 


0.49 


2.72 


65 Ge(p, 7 ) 66 As 


65 Zn 




3.74 


66 Ge(p, 7 ) 67 As 


66 Ge 


0.49 




67 As(p, 7 ) 68 Se 


^ Gel 


0.38 


3.20 


69 Se(p, 7 ) 70 Br 


69 Ge 


0.24 


3.37 




70 Ge 




0.44 


70 Se(p, 7 ) 71 Br 


70 Ge 


0.41 


2.66 


70 Br(p, 7 ) 71 Kr 


70 Ge 


0.33 




71 Br(p, 7 ) 72 Kr 


71 As 


0.24 


3.89 


73 Kr(p, 7 ) 74 Rb 


73 Se 


0.20 


3.64 




74 Se 




0.47 


74 Kr(p, 7 ) 75 Rb 


74 Se 


0.32 


2.91 


74 Rb(p, 7 ) 75 Sr 


74 Se 


0.41 
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Table 3 — Continued 



Reaction 


Isotope 


10 


0.1 


75 Rb(p, 7 ) 76 Sr 


75 Br 


0.21 


4.13 


76 Rb(p, 7 ) 77 Sr 


77 Kr 




0.48 


77 Sr(p, 7 ) 78 Y 


77 Kr 


0.19 


3.57 


78 Sr(p, 7 ) 79 Y 


78 Kr 


0.24 


3.18 


78 Y(p, 7 ) 79 Zr 


78 Kr 


0.48 




79 Y(p, 7 ) 80 Zr 


79 Kr 




4.99 


80 Y(p, 7 ) 81 Zr 


81 Rb 




0.43 


81 Zr(p, 7 ) 82 Nb 


81 Rb 


0.19 


3.35 


82 Zr(p, 7 ) 83 Nb 


82 Sr 


0.21 


3.21 




83 Sr 




0.46 


83 Nb(p, 7 ) 84 Mo 


83 Sr 




3.62 


84 Nb(p, 7 ) 85 Mo 


84g r 




4.18 


85 Mo(p, 7 ) 86 Tc 


85y 


0.19 


2.83 




90 Mo 




0.50 




91 Nb 




0.46 


86 Mo(p, 7 ) 87 Tc 


86 Zr 


0.17 


2.97 




87 Zr 




0.44 




88 Zr 




0.40 




89 Nb 




0.44 




90 Mo 




0.42 


87 Tc(p, 7 ) 88 Ru 


87 Zr 


0.21 






88 Zr 




0.43 


88 Tc(p, 7 ) 89 Ru 


89 Nb 




0.47 
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Table 4. Same as Table 1, for Model K04-B1. 



Reaction Isotope 10 0.1 



15 0(a, 7 ) 19 Ne 


15 N 


0.11 


6.14 


18 Nefa n , ) 21 Na 

J. i \ - yJ ) lift 


15 N 


0.28 






18p 


0.29 






21 Ne 


0.26 


2.03 




22 Na 


0.49 


6.21 


23 Al(p, 7 ) 24 Si 


22 Na 




4.45 


27 PfD V) 28 S 


27 A1 


0.25 


6.34 


31 Cl(p, 7 ) 32 Ar 


30 Si 


0.16 


2.57 




31p 


0.16 


2.72 




32 s 




0.29 


34 Ar(p, 7 ) 35 K 


34g 


0.44 


2.25 




36 Ar 




0.30 


39 Ca(p, 7 ) 40 Sc 


39 K 




9.70 




42 Ca 




2.11 


43 V(p, 7 ) 44 Cr 


42 Ca 


0.22 


7.20 


47 Mn(p, 7 ) 48 Fe 


46 Ti 


0.17 


2.38 




48 Cr 




0.26 


51 FefD V) 52 Co 


51 Cr 




6.19 




52 Fe 


0.35 


7.51 


52 Co(d ^ 53 Ni 


52 Fe 




2.24 


55 NifD V) 56 Cu 


55 Co 




2.83 


56 NifD V) 57 Cu 


56 Ni 


0.33 


7.74 




55 p„ 




a no 




56 Ni 




2.16 


57 Cu(p, 7 ) 58 Zn 


57 Ni 




2.33 


61 Ga(p, 7 ) 62 Ge 


60 Ni 


0.17 


2.93 




61 Cu 


0.16 


3.01 




62 Zn 




0.35 




64 Zn 




0.39 




65 Zn 




0.37 




66 Ge 




0.37 








0.38 




68 Ge 




0.33 
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Table 4 — Continued 



Reaction Isotope 10 0.1 





by Ge 




0.31 




70 Ge 




0.31 




71 As 




0.32 




72 Se 




0.30 




76 Kr 




0.29 


63 Ga(p, 7 ) 64 Ge 


63 Cu 




2.13 


65 Ge(p, 7 ) 66 As 


65 Zn 


0.25 


7.08 


65 As(p, 7 ) 66 Se 


68 Ge 


3.85 






69 Ge 


3.92 






70 Ge 


3.94 






71 As 


4.04 






72 Se 


4.73 


0.40 




76 Kr 


5.05 


0.32 


66 Ge(p, 7 ) 67 As 


66 Ge 




2.68 


66 As(p, 7 ) 67 Se 


66 Ge 


0.19 


2.66 


67 As(p, 7 ) 68 Se 


67 Ga 




3.15 


69 Se(p, 7 ) 70 Br 


69 Ge 




5.91 


70 Se(p, 7 ) 71 Br 


70 Ge 




2.96 


71 Br(p, 7 ) 72 Kr 


71 As 


0.42 


4.16 



Table 5. Same as Table 1, for Model K04-B2. 



F» portion 




in 


1 


12 C(a, 7 ) 16 


12 C 


0.19 


2.37 




24 Mg 




0.26 




28 Si 




0.29 


20 Ne(a, 7 ) 24 Mg 


20 Ne 




8.28 


24 Mg(a, 7 ) 28 Si 


24 Mg 


0.16 


3.77 




34 S 


0.48 


4.30 


24 Mg(a, p) 27 Al 


34 S 


3.99 


0.49 


28 Si(a, 7 ) 32 S 


28 Si 


0.31 


2.63 




34 S 




2.56 


32 S(a, 7 ) 36 Ar 


32 S 


0.31 






33 S 


0.48 






35 C1 


0.24 






36 Ar 


2.43 


0.35 




39 K 


2.20 




32 S(a, p) 35 Cl 


35 C1 


2.22 


0.30 


33 S(p, 7 ) 34 C1 


37 Ar 


0.43 




33 S(a, 7 ) 37 Ar 


33 S 


0.40 






37 Ar 


4.12 




34 S(p, 7 ) 35 C1 


34 S 


0.35 


2.02 


34 S(a, 7 ) 38 Ar 


34 S 


0.41 






38 Ar 


2.18 




35 Cl(p, 7 ) 36 Ar 


35 C1 


0.34 




35 Cl(a, p) 38 Ar 


38 Ar 


4.47 


0.45 


36 Ar(a, p) 39 K 


39 K 


5.47 


0.23 


39 K(p, 7 ) 40 Ca 


39 K 


0.45 


2.07 


40 Ca(a, p) 43 Sc 


44 Ti 




2.78 


41 Sc(p, 7 ) 42 Ti 


44 Ti 




0.50 


43 Sc(p, 7 ) 44 Ti 


44 Ti 


2.78 




51 Mn(p, 7 ) 52 Fe 


51 Cr 




3.66 


55 Co(p, 7 ) 56 Ni 


55 Co 




2.89 


60 Cu(p, 7 ) 61 Zn 


61 Cu 


4.92 


0.15 


61 Ga(p, 7 ) 62 Ge 


60 Ni 




2.54 




61 Cu 




2.12 
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Table 5 — Continued 



Reaction Isotope 10 0.1 





63 Cu 




2.56 


63 Ga(p, 7 ) 64 Ge 


63 Cu 




2.73 


64 Ga(p, 7 ) 65 Ge 


65 Zn 


5.48 


0.26 


65 As(p, 7 ) 66 Se 


36 Ar 




0.49 




39 K 




0.38 




61 Cu 




0.22 




64 Zn 


0.31 


6.16 




65 Zn 


0.38 


2.26 




66 Ge 


0.34 


4.24 




67 Ga 




11.28 




71 As 




0.33 




73 As 




0.42 




73 Se 




0.44 




84g r 




0.28 


/~""7 ---- / \ r- o > 

67 Ge(p, 7) 68 As 


67 Ga 


3.47 




67 As(p, 7 ) 68 Se 


66 Ge 




2.09 




67 Ga 




2.49 


68 As(p, 7 ) 69 Se 


69 Ge 


2.46 


0.42 


69 As(p, 7 ) 70 Se 


70 Ge 


2.36 




71 Se(p, 7 ) 72 Br 


71 As 


8.13 


0.13 


71 Br(p, 7 ) 72 Kr 


70 Ge 


0.44 


3.95 


72 Br(p, 7 ) 73 Kr 


73 As 


2.75 


0.37 




73 Se 


2.71 


0.37 


74 Kr(p, 7 ) 75 Rb 


75 Br 


0.28 




75 Rb(p, 7 ) 76 Sr 


74 Se 




2.77 




75 Se 




3.79 




75 Br 


0.24 


3.77 


76 Rb(p, 7 ) 77 Sr 


77 Br 


2.52 


0.46 




77 Kr 


2.49 


0.46 


78 Sr(p, 7 ) 79 Y 


78 Kr 




3.16 


79 Y(p, 7 ) 80 Zr 


79 Kr 


0.12 


5.75 


82 Zr(p, 7 ) 83 Nb 


82 Sr 




3.49 


83 Nb(p, 7 ) 84 Mo 


83 Sr 




2.50 



-45 - 



Table 5 — Continued 



Reaction 


Isotope 


10 


0.1 


84 Zr(p, 7 ) 85 Nb 


84 Sr 


2.35 


0.17 




85g r 


0.16 






85y 


0.15 




84 Nb(p, 7 ) 85 Mo 


83 Sr 




3.33 


86 Mo(p, 7 ) 87 Tc 


86y 




2.73 




86 Zr 




2.70 


87 Mo(p, 7 ) 88 Tc 


87y 


2.13 






87 Zr 


2.08 




88 Mo(p, 7 ) 89 Tc 


89 Zr 


0.19 






89 Nb 


0.19 




92 Ru(p, 7 ) 93 Rh 


92 Mo 


2.46 






93 Mn 

1V1U 


n 97 






93 Tc 


0.26 




92 Rh(p, 7 ) 93 Pd 


91 Nb 


2.14 




93 Rh(p, 7 ) 94 Pd 


92 Mo 




2.45 




93 Mo 




2.66 




93 Tc 




2.62 


96 Ag(p, 7 ) 97 Cd 


95 Tc 




2.06 




95 Ru 




2.03 




96 Ru 


0.31 


4.80 


"Cd(p, 7 ) 100 In 


"Rh 




3.19 


100 In(p, 7 ) 101 Sn 


ioo Rh 




3.12 




lOOpj 


0.31 


3.10 


101 Cd(p, 7 ) 102 In 


102p d 


0.45 
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Table 6. Same as Table 1, for Model K04-B3. 



Reaction Isotope 10 0.1 



15 0(a, 7) 19 Ne 


15 N 




8.17 


18 Ne(a, p) 21 Na 


15 N 


0.43 






18p 


0.44 




31 Cl(p, 7) 32 Ar 


30 Si 


0.45 




43 V(d Tl 44 Cr 
v IP) 1 j K - jL 


42 Ca 




4.68 




44 Ti 




0.48 


47 Mn(p 7) 48 Fe 


46 T j 


0.32 




56 Ni(p, 7) 57 Cu 


56 Ni 




7.93 


56 Cu(n ^ 57 Zn 


55 Co 




2.63 


57 Cu(d V) 58 Zn 


57 Ni 




2.79 


58 Cu(d t"1 59 Zii 


58 Ni 




2.08 


59 Cu(d V) 60 Zn 


59 Ni 




2.05 


61 ZnfD -/) 62 Ga 


61 Cu 




2.39 


61 GafD T) 62 Ge 


60 Ni 


0.17 


4.17 




61 Cu 




4.33 


63Q a f D ^)64q 


63 Cu 




3.31 




65 Zn 




8.88 




64 Zn 


0.34 


3.41 




65 Zn 


0.36 


2.81 




66 Ge 


0.43 


2.66 




67 Ga 


0.40 


3.12 




81 Rb 




0.48 




or 




n /la 




86 Zr 




0.50 




87 Zr 




0.49 




88 Zr 




0.48 




89 Nb 




0.48 




90 Mo 




0.47 




91 Nb 




0.46 




92 Mo 




0.46 




93 Tc 




0.45 




94 Mo 




0.44 




95 Ru 




0.44 
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Table 6 — Continued 



Reaction Isotope 10 0.1 





96 Ru 




0.43 




97 Ru 




0.41 




98 Ru 




0.39 




"Rh 




0.39 




ioo pd 




0.39 




ioi pd 




0.38 




102p d 




0.38 




103 Ag 




0.39 




104 Ag 




0.42 


66 Ge(p, 7 ) 67 As 


66 Ge 




4.81 


66 As(p, 7 ) 67 Se 


66 Ge 


0.18 




67 As(p, 7 ) 68 Se 


67 Ga 


0.45 


5.73 


69 Se(p, 7 ) 70 Br 


69 Ge 


0.10 


8.19 


70 Se(p, 7 ) 71 Br 


70 Ge 




5.24 


70 Br(p, 7 ) 71 Kr 


70 Ge 


0.28 




71 Br(p, 7 ) 72 Kr 


71 As 


0.20 


7.01 


73 Kr(p, 7 ) 74 Rb 


73 Se 


0.10 


7.53 


74 Kr p, 7 75 Rb 

\ L 7 / / 


74 Se 


0.30 


6.30 


74 Rb(p, 7 ) 75 Sr 


74 Se 


0.40 




75 Rb(p, 7 ) 76 Sr 


75 Br 


0.15 


6.45 


77 Sr(p, 7 ) 78 Y 


77 Kr 


0.11 


5.58 


78 Sr(p, 7 ) 79 Y 


78 Kr 


0.17 


6.15 


79 Y(p, 7 ) 80 Zr 


79 Kr 


0.12 


5.73 


81 Zr(p, 7 ) 82 Nb 


81 Rb 


0.15 


3.20 


82 Zr(p, 7 ) 83 Nb 


82 Sr 


0.17 


2.94 


83 Nb(p, 7 ) 84 Mo 


83 Sr 


0.29 




84 Nb(p, 7 ) 85 Mo 


84 Sr 


0.20 


2.51 


85 Nb(p, 7 ) 86 Mo 


85y 


0.44 


2.01 


85 Mo(p, 7 ) 86 Tc 


85y 


0.48 




86 Mo(p, 7 ) 87 Tc 


86 Zr 


0.23 


2.56 


87 Mo(p, 7 ) 88 Tc 


87 Zr 


0.32 




88 Tc(p, 7 ) 89 Ru 


88 Zr 


0.39 




89 Tc(p, 7 ) 90 Ru 


89 Nb 


0.25 


2.16 
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Table 6 — Continued 



Reaction 


Isotope 


10 


0.1 


90 Ru(p, 7 ) 91 Rh 


90 Mo 


0.33 




91 Ru(p, 7 ) 92 Rh 


91 Nb 


0.36 




92 Ru(p, 7 ) 93 Rh 


92 Mo 


0.32 




93 Rh(p, 7 ) 94 Pd 


93 Tc 


0.24 




94 Rh(p, 7 ) 95 Pd 


94 Mo 


0.45 




94 Pd(p, 7 ) 95 Ag 


94 Mo 


0.43 




95 Pd(p, 7 ) 96 Ag 


95 Ru 


0.30 




96 Ag(p, 7 ) 97 Cd 


96 Ru 


0.45 






97 Ru 




0.41 




98 Ru 




0.38 




"Rh 




0.43 




lOOpj 




0.48 


97 Ag(p, 7 ) 98 Cd 


97 Ru 


0.26 


2.04 


98 Ag(p, 7 )"Cd 


98 Ru 


0.45 




98 Cd(p, 7 )"In 


98 Ru 


0.44 




"Cd(p, 7 ) 100 In 


99 Rh 


0.33 




100 Cd(p, 7 ) 101 In 


ioo pd 


0.44 




101 In(p, 7 ) 102 Sn 


ioi pd 


0.34 




102 In(p, 7 ) 103 Sn 


102 pd 


0.38 






103 Ag 


2.05 


0.39 




104 Ag 




0.46 


103 In(p, 7 ) 104 Sn 


103 Ag 


0.50 






104 Ag 


2.81 


0.25 




105 Ag 


2.08 




104 In(p, 7 ) 105 Sn 


105 Ag 


4.13 
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Table 7. Same as Table 1, for Model K04-B4. 



Reaction Isotope 10 0.1 



12 C(a, 7 ) 16 


12 C 


0.24 


2.53 




16Q 


2.38 






20 Ne 


2.30 


0.26 




24 Mg 

b 




0.29 




28 Si 




0.41 


16 0(a V) 20 Ne 


16Q 




7.65 


20 Ne(a, 7) 24 Mg 


20 Ne 




7.87 


24 Me(a V) 28 Si 


24 Me 


0.20 


3.48 




30 Si 


0.44 


5.51 


24 M K (a, p) 27 Al 


30 Si 


5.84 


0.44 


28 Si(a p) 31 P 


31p 


3.39 


0.44 


28 Si(a o^ 32 S 


28 Si 


0.30 


3.04 




31p 


0.42 


4.07 


30 Si(a o^ 34 S 


30 Si 


0.39 


5.84 


30 P(p, 7 ) 31 S 


33 S 


0.50 






36 C1 


0.46 






30 Si 




23.47 




31p 




2.12 




36 C1 




0.37 


30 S(a, p) 33 Cl 


30 Si 




5.84 




31p 




2.88 




33 s 


0.04 


6.20 




36 pi 




ZL 1 Q 




37 Ar 


0.42 


2.56 


31 P(a, p) 34 S 


31p 




4.63 


32 S(a, 7 ) 36 Ar 


32 s 


0.36 






36 Ar 


2.35 


0.48 


33 S(a, p) 36 Cl 


36 C1 


5.59 




33 S(a, 7 ) 37 Ar 


33 S 


0.05 




34 S(a, 7 ) 38 Ar 


34 S 


0.23 




34 Cl(p, 7 ) 35 Ar 


34 S 


0.39 




34 Cl(a, p) 37 Ar 


34 S 


0.26 






37 Ar 


2.72 


0.44 
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Table 7 — Continued 



rtedCLion 


isoiope 


1 n 


n i 

U. 1 


34 Ar(a, p) 37 K 


34 s 


0.42 




35 Cl(a, p) 38 Ar 


35Q 


0.36 




37 Ar(a, 7 ) 41 Ca 


41 Ca 


3.93 




38 K(a, p) 41 Ca 


41 Ca 


6.45 


0.43 


39 K(p, 7 ) 40 Ca 


40 Ca 


2.28 




40 Ca(a, p) 43 Sc 


43 Ca 


0.40 






43 Sc 


0.40 




42 Ca(p, 7) 43 Sc 


42 Ca 


0.29 




43 Sc(p, 7 ) 44 Ti 


44 Ti 


2.99 




46 Ti(a, p) 49 V 


49y 


2.12 




48 Cr(p, 7 ) 49 Mn 


49y 




3.27 


49 Mn(p, 7) 50 Fe 


49y 




2.12 


58 Cu(p, 7) 59 Zn 


58 Ni 


6.90 


0.38 


60 Cu(p, 7 ) 61 Zn 


61 Cu 


3.59 


0.39 


61 Zn(p, 7 ) 62 Ga 


61 Cu 




2.41 


62 Zn(p, 7 ) 63 Ga 

\Jr 7 i / 


62 Zn 


4.80 


0.32 




63 Cu 


0.33 




64 Ga(p, 7 ) 65 Ge 


65 Zn 


2.52 




65 Ge(p, 7 ) 66 As 


65 Zn 




3.12 


65 As(p, 7) 66 Se 


64 Zn 


0.43 






66 Ge 


0.48 






^ Get 


0.48 




66 Ge(p, 7 ) 67 As 


66 Ge 




0.24 




67 Ga 




3.62 


70 Se(p, 7 ) 71 Br 


71 As 




2.65 


71 Se(p, 7 ) 72 Br 


71 As 


8.95 




73 Kr(p, 7 ) 74 Rb 


73 Se 




2.16 


74 Kr(p, 7 ) 75 Rb 


74 Se 




0.33 




75 Br 




4.28 


84 Zr(p, 7 ) 85 Nb 


85y 


0.44 




87 Mo(p, 7 ) 88 Tc 


88 & 


0.47 
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Table 8. Same as Table 1, for Model K04-B5. 



T-? DQpf in n 
rvcdCLlUll 




i n 


n 1 

U. 1 


12 C(a, 7 ) 16 


12 C 


0.26 






16Q 


2.54 






20 Ne 


2.44 


0.18 




24 Mg 


2.13 


0.21 




28 Si 




0.46 


16 0(a, 7 ) 20 Ne 


16Q 




7.66 


20 Ne(a, 7) 24 Mg 


20 Ne 




7.87 


24 Mg(a, 7 ) 28 Si 


24 Mg 


0.20 


3.38 


25 Si(a, p) 28 P 


30 Si 


2.12 


0.39 




32 S 


0.39 






36 Ar 


0.40 






37 Ar 


0.15 


2.54 




47 Ti 




2.10 




55 Co 




2.06 




56 Ni 


0.37 


2.10 




59 Ni 


0.36 


2.75 




65 Zn 


2.21 






66 Ge 


2.28 






68 Ge 


2.03 






69 Ge 


2.30 






70 Ge 


3.48 






71 As 


3.35 


0.35 




72 Se 


2.88 


0.48 




73 Se 


4.12 






76 Kr 


4.02 


0.36 




80 Sr 


5.57 




28 Si(p, 7 ) 29 P 


29 Si 


2.88 




28 Si(a, 7 ) 32 S 


28 Si 


0.30 


2.33 




32 S 




0.34 


29 Si(p, 7 ) 3 °P 


29 Si 


0.46 




29 Si(a, 7 ) 33 S 


29 Si 




3.62 


29 P(a, p) 32 S 


29 Si 




4.17 


29 S(a, p) 32 Cl 


32 S 


0.42 
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Table 8 — Continued 



Reaction Isotope 10 0.1 



30 P(p, 7 ) 31 s 
30 P(«, p) 33 S 
30 S(a, p) 33 Cl 



36 Ar 


0.41 




37 Ar 


0.18 




56 Ni 


0.38 




59 Ni 


0.34 




66 Ge 


2.02 




69 Ge 


2.02 




70 Ge 


2.84 




71 As 


2.81 




72 Se 


2.43 




73 Se 


3.39 




76 Kr 


3.30 




80 Sr 


4.51 




30 Si 


0.32 




31p 


2.11 


0.35 


30 Si 


0.03 


2.87 


33 s 




0.28 


32 s 


0.09 




33 s 


0.20 




34 s 


0.39 




35 CI 


0.26 




36 Ar 


0.17 


2.20 


37 Ar 


0.05 


3.72 


38 Ar 


0.36 


2.05 


39 K 


0.36 




40 Ca 


0.35 




42 Ca 


0.44 


2.18 


43 Sc 


0.41 




44 Ti 


0.33 




46 Ti 


0.38 


2.65 


47 Ti 


0.37 


2.63 


48 Ti 


0.37 


2.65 


50 Ti 


0.35 


2.65 


51 Ti 


0.36 


2.48 
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Table 8 — Continued 



Reaction 


IsotoDe 


10 


0.1 




54p e 


0.34 


2.67 




55 Co 


0.33 


2.63 




56 Ni 


0.23 






59 Ni 


0.25 


3.66 




61 Cu 


2.15 






62 Zn 


2.43 






65 Zn 


4.98 






66 Ge 


5.83 






67 Ga 


2.58 






68 Ge 


3.01 






69 Ge 


5.59 






70 Ge 


11.80 






71 As 


8.17 






72 Se 


5.36 






73 Se 


14.32 






76 Kr 


9.40 


0.49 




80g r 


15.91 




31 P(a, p) 34 S 


31 p 


0.17 




32 S(a, 7 ) 36 Ar 


36 Ar 


3.60 




33 S(a, 7 ) 37 Ar 


37 Ar 


5.33 


0.36 


34 Cl(a, p) 37 Ar 


37 Ar 


3.19 




38 K(p, 7 ) 39 Ca 


38 Ar 


0.44 




40 Ca(a, p) 43 Sc 


44 Ti 


0.48 




42 Ca(p, 7 ) 43 Sc 


42 Ca 


0.35 




43 Sc(p, 7 ) 44 Ti 


44 Ti 


4.48 


0.29 


46 Ti(p, 7 ) 47 V 


47 Ti 


2.65 




46 V(p, 7 ) 47 Cr 


47 Ti 


2.41 




48 Cr(p, 7 ) 49 Mn 


49y 


0.28 




49 Mn(p, 7 ) 50 Fe 


49y 




2.30 


52 Fe(p, 7 ) 53 Co 


52p e 




2.52 


53 Fe(p, 7 ) 54 Co 


53 Mn 




2.25 


53 Co(p, 7 ) 54 Ni 


53 Mn 




3.06 


56 Ni(a, p) 59 Cu 


35Q 


2.08 
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Table 8 — Continued 



Reaction Isotope 10 0.1 

36 Ar 2.74 ... 

37 Ar 3.74 ... 

46 Ti 2.18 ... 

47 Ti 2.70 ... 

48 Cr 2.30 ... 

50 Cr 2.24 ... 

51 Cr 2.18 ... 

54 Fe 2.29 ... 

55 Co 2.31 ... 



56 Ni 2.48 ... 

59 Ni 4.43 ... 

65 Zn 0.49 ... 

69 Ge 0.45 ... 

71 As 0.32 ... 

72 Se 0.42 ... 

76 Kr 0.31 ... 

57 Cu(p, 7 ) 58 Zn 56 Ni ... 2.85 

59 Ni 2.04 ... 

59 Cu(p, 7 ) 60 Zn 35 C1 ... 2.10 

36 Ar ... 2.79 

37 Ar ... 3.83 

46 Ti ... 2.20 

47 Ti ... 2.77 

48 Cr ... 2.35 

50 Cr ... 2.29 

51 Cr ... 2.21 

54 Fe ... 2.32 

55 Co ... 2.35 

56 Ni ... 2.53 

59 Ni ... 8.36 

65 Zn ... 0.48 

69 Ge ... 0.45 

71 As ... 0.30 
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Table 8 — Continued 



Reaction 


Isotope 


10 


0.1 




72 Se 




0.42 




76 Kr 




0.30 


61 Zn(p, 7 ) 62 Ga 


61 Cu 


0.25 


3.37 




62 Zn 




0.44 


61 Ga(p, 7 ) 62 Ge 


64 Zn 




0.47 




67 Ga 




0.49 




68 Ge 




0.46 




69 Ge 




0.47 


62 Ga(p, 7 ) 63 Ge 


62 Zn 


0.20 


2.75 


63 Ga(p, 7 ) 64 Ge 


63 Cu 




2.36 


64 Ga(p, 7 ) 65 Ge 


65 Zn 




0.50 


65 Ge(p, 7 ) 66 As 


65 Zn 


0.21 


3.68 


65 As(p, 7 ) 66 Se 


68 Ge 


2.11 






69 Ge 


2.10 






70 Gc 


2.14 






71 As 


2.15 






72 Se 


2.69 






oe 


z.uo 






76 Kr 


2.89 


0.47 




80 Sr 


2.97 




66 As(p, 7 ) 67 Se 


66 Ge 


0.29 


2.13 


67 As(p, 7 ) 68 Se 


67 Ga 


0.48 


2.78 


69 Se(p, 7 ) 70 Br 


69 Ge 




2.27 


70 Se(p, 7 ) 71 Br 


70 Gc 


2.91 






71 As 


0.38 




71 Br(p, 7 ) 72 Kr 


70 Ge 




2.36 




71 As 


0.30 


3.81 


73 Kr(p, 7 ) 74 Rb 


73 Se 




3.23 
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Table 9. Same as Table 1, for Model K04-B6. 



Reaction Isotope 10 0.1 



15 0(a, 7) 19 Ne 


15 N 




9.47 


18 Ne(a, p) 21 Na 


15 N 


0.42 






18p 


0.41 






22 Na 




3.14 


23 AKd V) 24 Si 


22 Na 




5.61 


27 P(p 7 )28g 


27 A1 




4.73 




28 Si 




0.48 


39 Ca(p, 7) 40 Sc 


39 K 




8.52 


43 V(p, 7) 44 Cr 


42 Ca 




2.65 


52 FefD 7) 53 Co 


52 Fe 




8.62 


52 Co(d Tl 53 Ni 


52 Fe 




2.20 


53 Co(d ^ 54 Ni 


53 Mn 




2.07 




56 Ni 


0.25 


8.48 


56 Cu(d V) 57 Zn 


55 Co 


0.30 




57 CWd V) 58 Zn 

VUI [_/ • I / £-111 


57 Ni 


0.39 


2.86 


58 Cu(d ^ 59 Zn 


58 Ni 




2.91 


59 Cu(d -/) 60 Zn 


59 Ni 




2.57 


61 ZnfD V) 62 Ga 


61 Cu 




3.74 


61 GafD V) 62 Ge 


60 Ni 


0.16 


5.05 




61 Cu 


0.15 


5.71 




63 Cu 




3.22 




69 Ge 




0.49 




73 c„ 

oe 




U.4o 




74 Se 




0.43 




75 Br 




0.43 




76 Kr 




0.47 




77 Kr 




0.40 




78 Kr 




0.40 




79 Kr 




0.40 




80 Sr 




0.44 




81 Rb 




0.39 




82 Sr 




0.40 




83g r 




0.40 
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Table 9 — Continued 



Reaction Isotope 10 0.1 





84g r 




0.41 




85y 




0.41 




86 Zr 




0.41 




87 Zr 




0.40 




88 Zr 




0.39 




89 Nb 




0.39 




90 Mo 




0.39 




91 Nb 




0.39 




92 Mo 




0.39 




93 Tc 




0.40 




94 Mo 




0.41 


63 Ga(p, 7 ) 64 Ge 


63 Cu 




4.97 


65 Ge(p, 7 ) 66 As 


65 Zn 


0.11 


9.19 


65 As(p, 7) 66 Se 


88 Zr 


2.01 






89 Nb 


2.06 






90 Mo 


2.12 






91 Nb 


2.19 






92 Mo 


2.25 






93 Tc 


2.27 






94 Mo 


2.30 




66 Ge(p, 7 ) 67 As 


66 Ge 


0.39 


6.27 


66 As(p, 7 ) 67 Se 


66 Ge 


0.23 




67 As(p, 7) 68 Se 




0.28 


7.20 


69 Se(p, 7 ) 70 Br 


69 Ge 


0.10 


7.68 


70 Se(p, 7 ) 71 Br 


70 Ge 


0.39 


6.29 


70 Br(p, 7 ) 71 Kr 


70 Ge 


0.39 




71 Br(p, 7 ) 72 Kr 


71 As 


0.14 


6.88 


73 Kr(p, 7 ) 74 Rb 


73 Se 


0.10 


6.83 


74 Kr(p, 7 ) 75 Rb 


74 Se 


0.22 


6.80 


75 Rb(p, 7 ) 76 Sr 


75 Br 


0.13 


5.73 


77 Sr(p, 7 ) 78 Y 


77 Kr 


0.13 


4.46 


78 Sr(p, 7 ) 79 Y 


78 Kr 


0.14 


5.35 


79 Y(p, 7 ) 80 Zr 


79 Kr 


0.12 


4.53 
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Table 9 — Continued 



Reaction Isotope 10 0.1 



80 Y(p, 7 ) 81 Zr 


81 Rb 




0.48 


81 Zr(p, 7 ) 82 Nb 


81 Rb 


0.19 


2.26 


82 Zr(p, 7 ) 83 Nb 


82 Sr 


0.21 


2.12 




87 Zr 




0.49 




88 Zr 




0.48 




89 Nb 




0.48 




90 Mo 




0.47 




91 Nb 




0.47 




92 Mo 




0.46 




93 Tc 




0.46 




94 Mo 




0.45 


83 Zr(p, 7 ) 84 Nb 


83 Sr 


0.49 




83 Nb(p, 7 ) 84 Mo 


83 Sr 


0.36 




84 Nb(p, 7 ) 85 Mo 


84 Sr 


0.29 






92 Mo 




0.48 




93 Tc 




0.44 




94 Mo 




0.41 


85 Nb(p, 7 ) 86 Mo 


85y 


0.36 




86 Mo(p, 7 ) 87 Tc 


86 Zr 


0.37 






91 Nb 




0.47 




92 Mo 




0.43 




93 Tc 




0.40 




94 Mo 




0.37 


87 Mo(p, 7 ) 88 Tc 


87 Zr 


0.34 






91 Nb 




0.48 




92 Mo 




0.45 




93 Tc 




0.43 




94 Mo 




0.40 


88 Mo(p, 7 ) 89 Tc 


88 Zr 


0.37 






89 Nb 




0.50 


89 Tc(p, 7 ) 90 Ru 


89 Nb 


0.31 




91 Tc(p, 7 ) 92 Ru 


91 Nb 


0.42 






92 Mo 




0.46 
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Table 9 — Continued 



Reaction 


Isotope 


10 


0.1 


91 Ru(p, 7 ) 92 Rh 


94 Mo 




0.49 


92 Ru(p, 7 ) 93 Rh 


92 Mo 


0.30 






93 Tc 




0.27 




94 Mo 




0.32 


93 Rh(p, 7 ) 94 Pd 


93 Tc 


0.38 




94 Rh(p, 7 ) 95 Pd 


94 Mo 


0.46 
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Table 10. Same as Table 1, for Model K04-B7. 



Pf\ ft inn 


L yj vj c 


m 


1 

W . J- 


15 0(a, 7 ) 19 Ne 


15 N 




7.89 


18 Ne(a, p) 21 Na 


18p 


0.48 




31 Cl(p, 7 ) 32 Ar 


30 Si 




2.17 


34 Ar(p, 7 ) 35 K 


34 s 


0.37 


2.14 


47 Mn(p, 7 ) 48 Fe 


46 Ti 




3.23 


57 Cu(p, 7 ) 58 Zn 


57 Ni 




2.12 


61 Ga(p, 7 ) 62 Ge 




0.13 


5.21 


65 As(p, 7 ) 66 Se 


64 Zn 


0.38 




67 As(p, 7 ) 68 Se 


67 Ga 




3.13 


69 Se(p, 7 ) 70 Br 


69 Ge 


0.10 


9.06 


69 Br(p, 7 ) 70 Kr 


68 Ge 


0.18 






69 Ge 


0.19 






70 Ge 


0.19 






71 As 


0.18 




70 Se(p, 7 ) 71 Br 


70 Ge 




3.39 


70 Br(p, 7 ) 71 Kr 


70 Ge 


0.17 


2.06 


71 Br(p, 7 ) 72 Kr 


71 As 


0.48 


5.63 


73 Kr(p, 7 ) 74 Rb 


73 Se 


0.10 


8.80 


74 Kr(p, 7 ) 75 Rb 


74 Se 




5.35 


74 Rb(p, 7 ) 75 Sr 


74 Se 


0.25 




75 Rb(p, 7 ) 76 Sr 


75 Br 


0.31 


6.68 


77 Sr(p, 7 ) 78 Y 


77 Kr 


0.11 


7.48 


78 Sr(p, 7 ) 79 Y 


78 Kr 


0.27 


6.46 


78 Y(p, 7 ) 79 Zr 


78 Kr 


0.46 




79 Y(p, 7 ) 80 Zr 


79 Kr 


0.13 


7.26 


81 Zr(p, 7 ) 82 Nb 


81 Rb 


0.12 


4.97 


82 Zr(p, 7 ) 83 Nb 


82 Sr 


0.14 


4.15 


83 Nb(p, 7 ) 84 Mo 


83 Sr 


0.17 


2.76 


84 Nb(p, 7 ) 85 Mo 


84 Sr 


0.14 


3.32 


85 Mo(p, 7 ) 86 Tc 


85y 


0.29 




86 Mo(p, 7 ) 87 Tc 


86 Zr 


0.17 


3.01 


87 Mo(p, 7 ) 88 Tc 


87 Zr 


0.43 




87 Tc(p, 7 ) 88 Ru 


87 Zr 


0.42 
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Table 10 — Continued 



rvedctioii 


isotope 


1 n 


n 1 

U.l 


88 Tc(p, 7) 89 Ru 


88 Zr 


0.30 


2.44 


89 Tc(p, 7) 90 Ru 


89 Nb 


0.28 


2.34 


90 Ru(p, 7 ) 91 Rh 


90 Mo 


0.25 


2.46 


91 Ru(p, 7 ) 92 Rh 


91 Nb 


0.32 




92 Ru(p, 7 ) 93 Rh 


92 Mo 


0.40 




92 Rh(p, 7 ) 93 Pd 


92 Mo 


0.45 




93 Rh(p, 7 ) 94 Pd 


93 Tc 


0.21 


2.41 


94 Pd(p, 7) 95 Ag 


94 Mo 


0.35 




95 Pd(p, 7) 96 Ag 


95 Ru 


0.26 


2.11 


96 Ag(p, 7 ) 97 Cd 


96 Ru 


0.39 


2.41 




97 Ru 




0.46 




98 Ru 




0.47 


97 Ag(p, 7 ) 98 Cd 


97 Ru 


0.25 


2.36 


97 Cd(p, 7) 98 In 


97 Ru 


0.47 




98 Ag(p, 7) 99 Cd 


98 Ru 


0.49 




98 Cd(p, 7 ) 99 In 


98 Ru 


0.32 




"Cdfp 7) 100 In 


99 Rh 


0.27 


2.10 


100 Cd(p, 7 ) 101 In 


ioo pd 


0.43 




100 In(p, 7 ) 101 Sn 


ioo pd 


0.41 




101 In(p, 7 ) 102 Sn 


ioi pd 


0.27 




102 In(p, 7 ) 103 Sn 


102 pd 


0.32 






103 Ag 




0.45 


103 In(p, 7 ) 104 Sn 


103 Ag 


0.45 






104 Ag 


2.01 


0.28 




105 Ag 




0.43 


104 In(p, 7 ) 105 Sn 


105 Ag 


3.03 


0.20 




106 Cd 




0.24 


105 Sn(p, 7 ) 106 Sb 


105 Ag 


0.43 






106 Cd 


2.16 


0.31 
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Table 11. Important /3-decay (and /3-delayed proton decay) rates, as found in our ten 

models. See text for details. 



Reaction 



Model K04 

64 Ge(/3+) 64 Ga 

68 Se(/3+) 68 As 

72 Kr(/?+) 72 Br 

76 Sr(/5+) 76 Rb 

80 Zr(/^+) 80 Y 

88 Ru(/3+) 88 Tc 

92 Pd(/? + ) 92 Rh 

"In(/?+)"Cd 

Model F08 

26 Si(/3+) 26m Al 

29 S(/?+) 29 P 

29 S(/3+ p) 28 Si 

30 S(/?+) 30 P 

33 Ar(/?+) 33 Cl 

33 Ar(/?+ p) 32 S 

34 Ar(/?+) 34 Cl 

37 Ca(/?+) 37 K 

37 Ca(/5+ p) 36 Ar 

38 Ca(/?+) 38 K 

60 Zn(/?+) 60 Cu 

68 Se(/3+) 68 As 

72 Kr(/?+) 72 Br 

76 Sr(/3+) 76 Rb 

80 Zr(/3+) 80 Y 

Model SOI ~ 

68 Se(/5+) 68 As 

72 Kr03+) 72 Br 

80 Zr(/?+) 80 Y 

88 Ru(/5+) 88 Tc 

92 Pd(/?+) 92 Rh 

101 Sn(/?+) 101 In 

102 Sn(/?+) 102 In 
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Table 11 — Continued 



Reaction 



103 Sn(/3+) 103 In 

104 Sn(/?+) 104 In 

105 Sn(/?+) 105 In 

106 Sn(/5+) 106 In 

106 SbQ3 + ) 106 Sn 

Model K04-B1 

25 Si(/?+) 25 Al 

25 Si(/?+ p) 24 Mg 

33 Ar(/5+) 33 Cl 

33 Ar(/?+ p) 32 S 

36 Ca(/?+ p) 35 Ar 

36 Ca(/5+) 36 K 

37 Ca(/3+) 37 K 

37 Ca(/?+ p) 36 Ar 

38 Ca(/?+) 38 K 

41 Ti(/3+) 41 Sc 

41 Ti(/5 + p) 40 Ca 

59 Zn(/5+) 59 Cu 

62 Ge(/3+) 62 Ga 

63 Ge(/3+) 63 Ga 

64 Ge(/?+) 64 Ga 

68 SeQ3+) 68 As 

Model K04-B2 

25 Si(/3+) 25 Al 

25 Si(/3+ p) 24 Mg 

68 Se(/3+) 68 As 

72 Kr(/3+) 72 Br 

76 Sr(/3+) 76 Rb 

80 Zr(/?+) 80 Y 

92 Pd(/5+) 92 Rh 

94 Pd(/5+) 94 Rh 

101 Sn(/3+) 101 In 

102 Sn(/?+) 102 In 
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Table 11 — Continued 



Reaction 



103 Sn(/3+) 103 In 

104 Sn(/3+) 104 In 

Model K04-B3 

68 Se(/5+) 68 As 

72 Kr(/?+) 72 Br 

76 Sr(/5+) 76 Rb 

80 Zr(/?+) 80 Y 

88 Ru(/5+) 88 Tc 

92 Pd(/3+) 92 Rh 

"In(/3+)"Cd 

Model K04-B4 

33 Ar(/3+) 33 Cl 

33 Ar(/?+ p) 32 S 

34 Ar(/3+) 34 Cl 

37 Ca(/3+) 37 K 

37 Ca(/?+ p) 36 Ar 

38 Ca(/?+) 38 K 

39 Ca(/?+) 39 K 

42 Ti(/5+) 42 Sc 

46 Cr(/3+) 46 V 

49 Fe(/3+) 49 Mn 

49 Fe(/5+ p) 48 Cr 

55 Ni(/3+) 55 Co 

58 Zn(/?+) 58 Cu 

68 Se(/3+) 68 As 

72 Kr(/?+) 72 Br 

76 Sr(/5 + ) 76 Rb 

80 Zr(/3+) 80 Y 

Model K04-B5 

21 Mg(/3+) 21 Na 

21 Mg(/5+ p) 20 Ne 

24 Si(/?+) 24 Al 

24 Si(/3+ p) 23 Mg 
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Table 11 — Continued 



Reaction 



25 Si(/3+) 25 Al 

25 Si(/?+ p) 24 Mg 

28 S(/3+) 28 P 

28 S(/3+ p) 27 Si 

29 S(/5+) 29 P 

29 S(/3+ p) 28 Si 

3O S03 + ) 3O P 

33 Ar(/?+) 33 Cl 

33 Ar(/3+ p) 32 S 

34 Ar(/3+) 34 Cl 

37 Ca(/?+) 37 K 

37 Ca(/?+ p) 36 Ar 

38 Ca(/?+) 38 K 

64 Ge(/3+) 64 Ga 

68 SeQ3+) 68 As 

Model K04-B6 

18 Ne(/3+) 18 F 

64 Ge(/?+) 64 Ga 

68 Se(/5+) 68 As 

72 Kr(/?+) 72 Br 

76 Sr(/5+) 76 Rb 

80 Zr(/3+) 80 Y 

Model K04-B7 

68 Se(/?+) 68 As 

72 Kr(/?+) 72 Br 

76 Sr(/3+) 76 Rb 

80 Zr(/?+) 80 Y 

88 Ru(/3+) 88 Tc 

92 Pd(/3+) 92 Rh 

101 Sn(/?+) 101 In 

102 Sn(/3+) 102 In 
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Table 12. List of reactions selected for the Q-value variation study, restricted to 
proton-capture reactions on A < 80 nuclei, with \Q\ < 1 MeV, and a Q-value uncertainty > 

50 keV (see Audi et al. 2003a, and Section 3). 



Reaction 


Q-value a (keV) 


Uncertainty 


21 Mg(p, 7 ) 22 A1 


17 


95 


24 Si(p, 7 ) 25 P 


-828 


197 


25 Si(p, 7 ) 26 P 


140 


196 


25 P(p, 7 ) 26 s 


191 


357 


26 P(p, 7 ) 27 S 


719 


281 


29 S(p, 7 ) 30 C1 


-314 


202 


30 S(p, 7 ) 31 Cl b 


294 


50 


30 Cl(p, 7 ) 31 Ar 


439 


284 


33 Ar(p, 7 ) 34 K 


-614 


298 


41 Ti(p, 7 ) 42 V 


-242 


220 


42 Ti(p, 7 ) 43 V 


192 


233 


45 Cr(p, 7 ) 46 Mn 


694 


515 


46 Cr(p, 7 ) 47 Mn 


78 


160 


49 Fe(p, 7 ) 50 Co 


-98 


224 


50 Fe(p, 7 ) 51 Co 


88 


161 


51 Fe(p, 7 ) 52 Co 


982 


67 


54 Ni(p, 7 ) 55 Cu 


-298 


302 


55 Ni(p, 7 ) 56 Cu 


555 


140 


58 Zn(p, 7 ) 59 Ga 


-888 


175 


59 Zn(p, 7 ) 60 Ga 


27 


118 


59 Ga(p, 7 ) 60 Ge 


936 


287 


60 Zn(p, 7 ) 61 Ga b 


192 


54 


63 Ge(p, 7 ) 64 As 


-100 


300 


64 Ge(p, 7 ) 65 As 


-80 


300 


67 Se(p, 7 ) 68 Br 


-560 


300 


68 Se(p, 7 ) 69 Br 


-450 


100 


72 Kr(p, 7 ) 73 Rb 


-600 


150 


75 Sr(p, 7 ) 76 Y 


-629 


549 


76 Sr(p, 7 ) 77 Y 


-50 


50 


76 Y(p, 7 ) 77 Zr c 


802 


400 


80 Zr(p, 7 ) 81 Nb 


-750 


150 



(keV) 
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a All Q- values and errors are from Audi et al. (2003a), except where indicated. They have 
been estimated in that work from systematic trends. 

b Reactions with experimental Q- values and errors in Audi et al. (2003a). 

c Q-value as used in the Hauser-Feshbach calculation for this rate in our network. No Q- 
value was estimated in Audi et al. (2003a) for this reaction. The associated uncertainty was 
assumed to be 400 keV. 



Table 13. Final abundance ratios, Xi/X ijStd , for Model K04 (Koike et al. 2004), resulting 
from variations on Q- values (see Table 12, and text). 



Reaction Isotope High-Q (Q+AQ) Low-Q (Q-AQ) 



26 P(p, 7 ) 27 S 


25 Mg 


0.35 




46 Cr(p, 7 ) 47 Mn 


46 Ti 


0.23 




55 Ni(p, 7 ) 56 Cu 


55 Co 




3.95 


60 Zn(p, 7 ) 61 Ga 


eo Ni 


0.47 




64 Ge(p, 7 ) 65 As 


64 Zn 


0.08 


4.13 




65 Zn 




3.04 




66 Ge 




3.05 




67 Ga 




3.65 




73 Se 




0.49 




74 Se 




0.50 




77 Kr 




0.42 




78 Kr 




0.42 








0.42 




80 Sr 




0.47 




81 Rb 




0.40 




82 Sr 




0.41 




83 Sr 




0.42 




84g r 




0.42 




85y 




0.42 




86 Zr 




0.40 




87 Zr 




0.39 




88 Zr 




0.38 




89 Nb 




0.38 




90 Mo 




0.37 




91 Nb 




0.36 




92 Mo 




0.36 




93 Tc 




0.36 




94 Mo 




0.35 




95 Ru 




0.34 




96 Ru 




0.33 




97 Ru 




0.31 




98 Ru 




0.30 




"Rh 




0.29 



Table 13 — Continued 



Reaction Isotope High-Q (Q+AQ) Low-Q (Q-AQ) 



10 °pd 


0.30 


101 pd 


0.29 


102 pd 


0.29 


103 Ag 


0.30 


104 Ag 


0.32 



Table 14. Same as Table 13, for Model F08. 



Reaction Isotope High-Q (Q+AQ) Low-Q (Q-AQ) 



45 Cr(p, 7 ) 46 Mn 45 Ti 



0.50 
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Table 15. Correlations between final yields and enhancement factors applied to reaction 
rates in the Monte Carlo approach, for Model K04. All correlations with r > 0.5 (see text) 
are included. Reactions that have an impact of less than a factor of 2 on the final yields of 
correlated isotopes (for enhancement factors ranging between roughly 0.1 and 10) are 
indicated by brackets in the last column, in which slopes of the different linear fits (final 
isotopic yield vs. rate enhancement factor, see e.g. figs. 8 and 9) are given. The error 
associated with the slopes is typically ~1 - 10%. 



Isotope 


Reaction 


Correlation coefficient, r 


olope 


15 N 


15 0(a, 7 ) 19 Ne 


-0.778 


-0.927 


30 Si 


31 Cl(p, 7 ) 32 Ar 


-0.751 


[-0.205] 


33 s 


18 Ne(a, p) 21 Na 


-0.601 


[-0.063] 


34 s 


34 Arfp V) 35 K 


-0.603 


[-0.148] 


37 Ar 


18 Ne(a, p) 21 Na 


-0.649 


[-0.061] 


38 Ar 


18 Ne(a, p) 21 Na 


-0.571 


[-0.062] 


44 Ti 


43 V(p, 7 ) 44 Cr 


0.648 


[0.134] 


46 Ti 


47 Mn(p, 7 ) 48 Fe 


-0.822 


-0.292 


50 Cr 


18 Ne(a, p) 21 Na 


-0.582 


[-0.061] 


53 Mn 


53 Co(p, 7) 54 Ni 


-0.628 


[-0.121] 

L J 


54 Fe 


18 Ne(a, p) 21 Na 


-0.507 


[-0.055] 


55 Co 


56 Cu(p, 7 ) 57 Zn 


-0.777 


-0.484 


56 Ni 


56 Ni(p, 7) 57 Cu 


-0.840 


-0.703 


57 Ni 


57 Cu(p, 7 ) 58 Zn 


-0.807 


-0.376 


58 Ni 


58 Cu(p, 7) 59 Zn 


-0.593 


[-0.176] 


59 Ni 


59 Cu(p, 7 ) 60 Zn 


-0.658 


[-0.172] 

L J 




61 Ga(p, 7) 62 Ge 


-0.910 


-0.688 


61 Cu 


61 Ga(p, 7 ) 62 Ge 


-0.893 


-0.701 


62 Zn 


62 Ga(p, 7 ) 63 Ge 


-0.645 


[-0.185] 




61 Ga(p, 7 ) 62 Ge 


0.620 


[0.178] 


63 Cu 


63 Ga(p, 7 ) 64 Ge 


-0.779 


[-0.327] 


64 Zn 


65 As(p, 7 ) 66 Se 


-0.960 


-0.517 


65 Zn 


65 Ge(p, 7 ) 66 As 


-0.848 


-0.886 


66 Ge 


66 As(p, 7 ) 67 Se 


-0.585 


-0.547 




66 Ge(p, 7 ) 67 As 


-0.566 


-0.528 


67 Ga 


67 As(p, 7 ) 68 Se 


-0.699 


-0.549 




65 As(p, 7 ) 66 Se 


-0.553 


-0.440 


68 Ge 


69 Br(p, 7 ) 70 Kr 


-0.697 


[-0.063] 


69 Ge 


69 Se(p, 7 ) 70 Br 


-0.972 


-0.961 


70 Ge 


70 Se(p, 7 ) 71 Br 


-0.754 


-0.492 


71 As 


71 Br(p, 7 ) 72 Kr 


-0.963 


-0.795 


72 Se 


65 As(p, 7 ) 66 Se 


0.748 


[0.079] 


73 Se 


73 Kr(p, 7 ) 74 Rb 


-0.969 


-0.943 
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Table 15 — Continued 



Isotope 


Reaction 


Correlation coefficient, r 


Slope 


74 Se 


74 Kr(r> M 75 Rb 


-D 89^ 


-0 699 


75 o 

J_)l 


rvui \jj j j oi 




n si 7 


76 K 




n 71 1 

U. 1 11 


If) 1 1 91 


77 Kr 


° r lP) 7J 1 


-U. c/UZr 


-yj.o^to 


XYl 


78c r( ' r . ^79^ 


-U. jOO 


D 7Q/L 


79 Kr 

IY1 


79 Yfr> sy) m 7,r 
1 IP; 7J ZjI 




-u. o^to 


Ol 


80 Yfn Vl 81 Zr 




\ 1 Q4l 




8i7 r f n ^ 82 Nb 
^HP) 7J iNU 






82 c r 

Ol 


82 Zrfn -vl 83 Nb 


n qi s 


n fu^ 

U.U^iO 


83 c 

Ol 


83 Nbfn o/^Mn 

IN Ul ^ J 1V1U 


n 7^n 


D 498 


84 c r 


IN D^p, 7 J IVIO 


n 8^9 
-u.ooz 


-u.ooo 


85v 


85atu/„ A86i\/f„ 
IN u^p, 7 J IVIO 


-u.ooo 


-U.O 1 1 


867; 
ZjI 


ivio^p, 7j ±c 


n 81 8 

U.O 1 o 


D ^1 Q 


87 Zr 


87 Mo(p, 7 ) 88 Tc 


-0.752 


-0.412 


88 & 


88 Tc(p, 7 ) 89 Ru 


-0.611 


-0.360 


89 Nb 


89 Tc(p, 7 ) 90 Ru 


-0.755 


-0.467 


90 Mo 


90 Ru(p, 7 ) 91 Rh 


-0.680 


-0.391 


91 Nb 


91 Ru(p, 7 ) 92 Rh 


-0.583 


-0.329 


92 Mo 


92 Ru(p, 7 ) 93 Rh 


-0.699 


-0.429 


93 Tc 


93 Rh(p, 7 ) 94 Pd 


-0.716 


-0.460 


94 Mo 


94 Rh(p, 7 ) 95 Pd 


-0.585 


-0.358 


95 Ru 


95 Pd(p, 7 ) 96 Ag 


-0.648 


-0.388 


96 Ru 


96 Ag(p, 7 ) 97 Cd 


-0.591 


-0.331 


97 Ru 


97 Ag(p, 7 ) 98 Cd 


-0.549 


-0.420 


104 Ag 


103 In(p, 7 ) 104 Sn 


0.597 


0.589 


105 Ag 


104 In(p, 7 ) 105 Sn 


0.550 


0.623 



Table 16. Isotopes most affected by nuclear uncertainties, according to our sensitivity 
studies. Results from the individual reaction-rate approach (Table 1) are compared with 
those from Monte Carlo simulations (Table 15) for Model K04 (see text for details). 



Isotope Reaction (Individual variations) Reaction (Monte Carlo) 



15 N 15 0(a, 7 ) 19 Ne 15 0(a, 7 ) 19 Ne 

18 Ne(a, p) 21 Na 
18 F 18 Ne(a, p) 21 Na 

30 Si 31 Cl(p, 7 ) 32 Ar 

33 S 

44 Ti 43 V(p, 7 ) 44 Cr 

46 Ti 47 Mn(p, 7 ) 48 Fe 47 Mn(p, 7 ) 48 Fe 

55 Co 56 Cu(p, 7 ) 57 Zn 56 Cu(p, 7 ) 57 Zn 

56 Ni 56 Ni(p, 7 ) 57 Cu 56 Ni(p, 7 ) 57 Cu 

57 Ni 57 Cu(p, 7 ) 58 Zn 57 Cu(p, 7 ) 58 Zn 

58 Ni 58 Cu(p, 7 ) 59 Zn 

59 Ni 59 Cu(p, 7 ) 60 Zn 

60 Ni 61 Ga(p, 7 ) 62 Ge 61 Ga(p, 7 ) 62 Ge 

61 Cu 61 Ga(p, 7 ) 62 Ge 61 Ga(p, 7 ) 62 Ge 

61 Zn(p, 7 ) 62 Ga 
63 Cu 63 Ga(p, 7 ) 64 Ge 

64 Zn 65 As(p, 7 ) 66 Se 65 As(p, 7 ) 66 Se 

65 Zn 65 Ge(p, 7 ) 66 As 65 Ge(p, 7 ) 66 As 

65 As(p, 7 ) 66 Se 

66 Ge 66 Ge(p, 7 ) 67 As 66 Ge(p, 7 ) 67 As 

65 As(p, 7 ) 66 Se 

66 As(p, 7 ) 67 Se 66 As(p, 7 ) 67 Se 

67 Ga 67 As(p, 7 ) 68 Se 67 As(p, 7 ) 68 Se 

65 As(p, 7 ) 66 Se 65 As(p, 7 ) 66 Se 

69 Ge 69 Se(p, 7 ) 70 Br 69 Se(p, 7 ) 70 Br 

70 Ge 70 Se(p, 7 ) 71 Br 70 Se(p, 7 ) 71 Br 

70 Br(p, 7 ) 71 Kr 

71 As 71 Br(p, 7 ) 72 Kr 71 Br(p, 7 ) 72 Kr 

73 Se 73 Kr(p, 7 ) 74 Rb 73 Kr(p, 7 ) 74 Rb 

74 Se 74 Kr(p, 7 ) 75 Rb 74 Kr(p, 7 ) 75 Rb 

74 Rb(p, 7 ) 75 Sr 

75 Br 75 Rb(p, 7 ) 76 Sr 75 Rb(p, 7 ) 76 Sr 

77 Kr 77 Sr(p, 7 ) 78 Y 77 Sr(p, 7 ) 78 Y 



Table 16 — Continued 



Isotope Reaction (Individual variations) Reaction (Monte Carlo) 



78 Kr 


78 Sr(p, 7 ) 79 Y 


78 Sr(p, 7 ) 79 Y 


79 Kr 


79 Y(p, 7 ) 80 Zr 


79 Y(p, 7 ) 80 Zr 


81 Rb 


81 Zr(p, 7 ) 82 Nb 


81 Zr(p, 7 ) 82 Nb 




65 As(p, 7 ) 66 Se 


_ 


82 Sr 


82 Zr(p, 7 ) 83 Nb 


82 Zr(p, 7 ) 83 Nb 


83 Sr 


83 Nb(p, 7 ) 84 Mo 


83 Nb(p, 7 ) 84 Mo 


84 Sr 


84 Nb(p, 7 ) 85 Mo 


84 Nb(p, 7 ) 85 Mo 


85y 


85 Mo(p, 7 ) 86 Tc 






85 Nb(p, 7 ) 86 Mo 


85 Nb(p, 7 ) 86 Mo 


86 Zr 


86 Mo(p, 7 ) 87 Tc 


86 Mo(p, 7 ) 87 Tc 


87 Zr 


87 Mo(p, 7 ) 88 Tc 


87 Mo(p, 7 ) 88 Tc 


88 Zr 


88 Tc(p, 7 ) 89 Ru 


88 Tc(p, 7 ) 89 Ru 




65 As(p, 7 ) 66 Se 




89 Nb 


89 Tc(p, 7 ) 90 Ru 


89 Tc(p, 7 ) 90 Ru 




65 As(p, 7 ) 66 Se 




90 Mo 


90 Ru(p, 7 ) 91 Rh 


90 Ru(p, 7 ) 91 Rh 

\± 7 1 / 




65 As(p, 7 ) 66 Se 




91 Nb 


91 Ru(p, 7 ) 92 Rh 


91 Ru(p, 7 ) 92 Rh 




65 As(p, 7 ) 66 Se 


_ 


92 Mo 


92 Ru(p, 7 ) 93 Rh 


92 Ru(p, 7 ) 93 Rh 




65 As(p, 7 ) 66 Se 


_ 


93 Tc 


93 Rh(p, 7 ) 94 Pd 


93 Rh(p, 7 ) 94 Pd 




65 As(p, 7 ) 66 Se 


_ 


94 Mo 


65 As(p, 7 ) 66 Se 






94 Pd(p, 7 ) 95 Ag 






94 Rh(p, 7 ) 95 Pd 


94 Rh(p, 7 ) 95 Pd 


95 Ru 


95 Pd(p, 7 ) 96 Ag 


95 Pd(p, 7 ) 96 Ag 




65 As(p, 7 ) 66 Se 




96 Ru 


65 As(p, 7 ) 66 Se 






96 Ag(p, 7 ) 97 Cd 


96 Ag(p, 7 ) 97 Cd 


97 Ru 


97 Ag(p, 7 ) 98 Cd 


97 Ag(p, 7 ) 98 Cd 




65 As(p, 7 ) 66 Se 






96 Ag(p, 7 ) 97 Cd 
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Table 16 — Continued 



isot ope 


. . . 

Reaction (Individual variations) 


rteacxion ^ivionce oario ) 


98 Ru 


96 Ag(p, 7 ) 97 Cd 






65 As(p, 7 ) 66 Se 






98 Cd(p, 7 )"ln 






98 Ag(p, 7 )"Cd 




"Rh 


"Cd(p, 7 ) 100 In 






65 As(p, 7 ) 66 Se 






96 Ag(p, 7 ) 97 Cd 


- 


ioo pd 


65 As(p, 7 ) 66 Se 






100 Cd(p, 7 ) 101 In 






96 Ag(p, 7 ) 97 Cd 




ioi pd 


101 In(p, 7 ) 102 Sn 






65 As(p, 7 ) 66 Se 




102p d 


102 In(p, 7 ) 103 Sn 






65 As(p, 7 ) 66 Se 




103 Ag 


65 As(p, 7 ) 66 Se 






102 In(p, 7 ) 103 Sn 




104 Ag 


103 In(p, 7 ) 104 Sn 


103 In(p, 7 ) 104 Sn 




65 As(p, 7 ) 66 Se 






102 In(p, 7 ) 103 Sn 




105 Ag 


104 In(p, 7 ) 105 Sn 


104 In(p, 7 ) 105 Sn 




103 In(p, 7 ) 104 Sn 
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Table 17. Same as Table 15, for Model F08. 







( nrypl q l on r*/"i£i in pianf t* 
vjUI I cldllUll UUtJlllCltJIlL j 1 


Slnnp 

QlUjJL, 


4 He 


12 C(a, 7 ) 16 


-0.679 


[-0.030] 


12 C 


12 C(a, 7 ) 16 


-0.585 


[-0.203] 




12 C(p, 7 ) 13 N 


-0.581 


[-0.203] 


16Q 


16 0(a, 7 ) 20 Ne 


-0.878 


-0.945 


20 Ne 


20 Ne(a, 7 ) 24 Mg 


-0.894 


-1.01 


24 Mg 


24 Mg(a, 7 ) 28 Si 


-0.829 


-0.667 


26 Mg 


26 °A\(a, p) 29 Si 


-0.686 


-0.631 




26 Si(a, p) 29 P 


-0.550 


-0.493 


26 9A1 


26 °A\(a, p) 29 Si 


-0.709 


-0.788 




26 Si(a, p) 29 P 


-0.540 


-0.584 


29 Si 


2 ^Al(a, p) 29 Si 


0.570 


[0.221] 


30 Si 


30 P(p, 7 ) 31 S 


-0.687 


[-0.141] 


31p 


30 P(p, 7 ) 31 S 


0.856 


0.466 


32 s 


28 Si(a, 7 ) 32 S 


0.822 


0.467 


33 s 


30 P(a, p) 33 S 


0.587 


[0.204] 


34 s 


34 Cl(p, 7 ) 35 Ar 


-0.590 


[-0.108] 


35 C1 


34 Cl(p, 7 ) 35 Ar 


0.661 


[0.283] 


37 Ar 


36 Ar(p, 7 ) 37 K 


-0.915 


-0.482 


38 Ar 


38 K(p, 7 ) 39 Ca 


-0.606 


[-0.084] 


39 K 


40 Sc(p, 7 ) 41 Ti 


-0.704 


[-0.142] 


40 Ca 


39 K(p, 7 ) 40 Ca 


0.594 


[0.191] 


43 Sc 


44 V(p, 7 ) 45 Cr 


-0.838 


[-0.249] 


44 Ti 


44 V(p, 7 ) 45 Cr 


0.540 


[0.190] 


45 T j 


44 Ti(p, 7 ) 45 V 


-0.782 


-0.593 


47 Ti 


47 Cr(p, 7 ) 48 Mn 


-0.745 


-0.417 


48 Cr 


48 Cr(p, 7 ) 49 Mn 


-0.646 


[-0.173] 


49y 


49 Mn(p, 7 ) 50 Fe 


-0.859 


-0.295 


50 Cr 


50 Mn(p, 7 ) 51 Fe 


-0.534 


[-0.082] 


51 Cr 


52 Co(p, 7 ) 53 Ni 


-0.769 


[-0.162] 


52 Fe 


52 Fe(p, 7 ) 53 Co 


-0.934 


-0.513 


53 Mn 


53 Co(p, 7 ) 54 Ni 


-0.918 


-0.675 


54p e 


54 Co(p, 7 ) 55 Ni 


-0.721 


[-0.253] 


55 Co 


55 Co(p, 7 ) 56 Ni 


-0.761 


[-0.183] 
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Table 17— Continued 



Isotope Reaction Correlation coefficient, r Slope 



56 Ni 


57 Cu(p, 7 ) 58 Zn 


-0.805 


[-0.2651 

L J 


57 Ni 


57 Ni(p, 7 ) 58 Cu 


-0.900 


-0.858 


58 Ni 


58 Cu(p, 7 ) 59 Zn 


-0.754 


-0.584 




57 Cu p, 7 58 Zn 


0.624 


0.478 


59 Ni 


57 Cu(p, 7 ) 58 Zn 


0.670 


0.328 




59 Cu(p, 7 ) 60 Zn 


-0.668 


-0.327 




57 Cu(p, 7 ) 58 Zn 


0.614 


[0.096] 


61 Cu 


61 Zn(p, 7 ) 62 Ga 

\i it/ 


-0.902 


-0.621 


62 Zn 


62 Ga(p, 7 ) 63 Ge 


-0.777 


-0.498 


63 Cu 


63 Ga(p, 7 ) 64 Ge 


-0.865 


-0.406 


64 Zn 


61 Ga(p, 7 ) 62 Ge 


0.938 


0.345 


65 Zn 


65 Ge(p, 7 ) 66 As 


-0.831 


-0.662 


66 Ge 


66 As(p, 7 ) 67 Se 


-0.621 


-0.462 


67 Ga 


67 As(p, 7 ) 68 Se 


-0.747 


-0.495 




61 Ga(p, 7 ) 62 Ge 


0.528 


[0.349] 


68 Ge 


61 Ga(p, 7 ) 62 Ge 


0.851 


[0.031] 


69 Ge 


69 Se(p, 7 ) 70 Br 


-0.945 


-0.599 


70 Ge 


70 Se(p, 7 ) 71 Br 


-0.676 


-0.427 




70 Br(p, 7 ) 71 Kr 


-0.518 


-0.330 


71 As 


71 Br(p, 7 ) 72 Kr 


-0.935 


-0.627 


72 Se 


72 Br(p, 7 ) 73 Kr 


-0.514 


[-0.005] 

L J 


73 Se 


73 Kr(p, 7 ) 74 Rb 


-0.936 


-0.654 


74 Se 


74 Kr(p, 7 ) 75 Rb 


-0.797 


-0.514 


75 Br 


75 Rb(p, 7 ) 76 Sr 


-0.938 


-0.660 


76 Kr 


75 Rb(p, 7 ) 76 Sr 


0.571 


[0.020] 


77 Kr 


77 Sr(p, 7 ) 78 Y 


-0.926 


-0.653 


78 Kr 


78 Sr(p, 7 ) 79 Y 


-0.862 


-0.573 


79 Kr 


79 Y(p, 7 ) 80 Zr 


-0.932 


-0.771 


81 Rb 


81 Zr(p, 7 ) 82 Nb 


-0.903 


-0.639 


82 Sr 


82 Zr(p, 7 ) 83 Nb 


-0.912 


-0.596 


83 Sr 


83 Nb(p, 7 ) 84 Mo 


-0.916 


-0.679 


84 Sr 


84 Nb(p, 7 ) 85 Mo 


-0.935 


-0.750 


85y 


85 Mo(p, 7 ) 86 Tc 


-0.901 


-0.569 



Table 17— Continued 



Isotope 


Reaction 


Correlation coefficient, r 


Slope 


86 Zr 


86 Mo(p, 7 ) 87 Tc 


-0.916 


-0.617 


87 Zr 


87 Tc(p, 7 ) 88 Ru 


-0.752 


-0.451 


88 Zr 


87 Tc(p, 7 ) 88 Ru 


0.526 


[0.289] 
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Table 18. Same as Table 16, for Model F08. 



Isotope Reaction (Individual variations) Reaction (Monte Carlo) 



12 C 


12 C(p, 7 ) 13 N 

\L 7 1 / 






l2 C(a, 7 ) 16 




16Q 


16 0(a, 7 ) 20 Ne 


16 0(a, 7 ) 20 Ne 




12 C(a, 7 ) 16 




20 Ne 


20 Ne(a, 7) 24 Mg 


20 Ne(a, 7) 24 Mg 




12 C(a, 7 ) 16 




24 Mg 


24 Mg(a T) 28 Si 


24 Mg(a 7) 28 Si 




12 C(a, 7 ) 16 




25 Mg 

o 


22 Mg(a, p) 25 Al 






24 Mg(a, 7 ) 28 Si 






25 Al(a, p) 28 Si 






24 Mg(p, 7) 25 A1 




o 


2% s Al( a , p ) 29 Si 


26 »Al(a, p) 29 Si 






26 Si(a, p) 29 P 




263 Al(a, p) 29 Si 


269 M(a, p) 29 Si 




26 Si(a p) 29 P 


26 Si(a p) 29 P 


27 A1 


26 Si(a p) 29 P 






27 Si(a, p) 30 P 






27 Al(a, p) 30 Si 






26 » Al(p, 7 ) 27 Si 

\1 7 1/ 






26 9M{a, p) 29 Si 






22 Mg(a, p) 25 Al 




28 Si 


28 Si(a, 7 ) 32 S 




31p 


30 P(p, 7 ) 31 S 


30 P(p, 7 ) 31 S 


32 s 


28 Si(a, 7 ) 32 S 


28 Si(a, 7 ) 32 S 


33 s 


30 P(a, p) 33 S 






29 Si(a, 7 ) 33 S 




35Q 


34 Cl(p, 7 ) 35 Ar 




37 Ar 


36 Ar(p, 7 ) 37 K 


36 Ar(p, 7 ) 37 K 


40 Ca 


39 K(p, 7 ) 40 Ca 




45 Ti 


44 Ti(p, 7 ) 45 V 


44 Ti(p, 7 ) 45 V 




44 V(p, 7 ) 45 Cr 






45 V(p, 7 ) 46 Cr 
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Table 18— Continued 



Isotope Reaction (Individual variations) Reaction (Monte Carlo) 



47 Ti 


47 Cr(p, 7) 48 Mn 


47 Cr(p, 7) 48 Mn 




46 V(p, 7) 47 Cr 




49 y 


49 Mn(p, 7) 50 Fe 


49 Mn(p, 7) 50 Fe 


52 Fe 


52 Fe(p, 7 ) 53 Co 


52 Fe(p, 7) 53 Co 


53 Mn 


53 Co(p, 7 ) 54 Ni 


53 Co(p, 7 ) 54 Ni 




cr o -, — , / \ cro 

52 Fe(p, 7) 53 Co 


— 


54 Fe 


54 Co(p, 7) 55 Ni 


— 


56 Ni 


57 Cu(p, 7 ) 58 Zn 




57 Ni 


57 Ni(p, 7 ) 58 Cu 


57 Ni(p, 7 ) 58 Cu 




57 Cu(p, 7) 58 Zn 




58 Ni 


58 Cu(p, 7 ) 59 Zn 


pro / \ rr n . . 

58 Cu(p, 7 ) 59 Zn 




57 Cu(p, 7 ) 58 Zn 


57 Cu(p, 7 ) 58 Zn 


59 Ni 


59 Cu(p, 7) 60 Zn 


rn , / \ f'c\ r—, 

59 Cu(p, 7) 60 Zn 




cr T ~- / \ cr o - . 

57 Cu(p, 7 ) 58 Zn 


cr >7 / \ cr o . 

57 Cu(p, 7 ) 58 Zn 


61 Cu 


61 Zn(p, 7) 62 Ga 


61 Zn(p, 7) 62 Ga 


62 Zn 


62 Ga(p, 7) 63 Ge 


62 Ga(p, 7 ) 63 Ge 




61 Zn(p, 7) 62 Ga 




63 Cu 


63 Ga(p, 7) 64 Ge 


63 Ga(p, 7) 64 Ge 


64 Zn 


61 Ga(p, 7) 62 Ge 


61 Ga(p, 7) 62 Ge 


65 Zn 


a cr / \ eta i 

65 Ge(p, 7) 66 As 


c cr • / \ act i 

65 Ge(p, 7) 66 As 




61 Ga(p, 7 ) 62 Ge 


— 


66 Ge 


61 Ga(p, 7 ) 62 Ge 
66 Ge(p, 7 ) 67 As 


— 






66 As(p, 7 ) 67 Se 


67 Ga 


67 As(p, 7 ) 68 Se 


67 As(p, 7 ) 68 Se 




61 Ga(p, 7 ) 62 Ge 




69 Ge 


69 Se(p, 7 ) 70 Br 


69 Se(p, 7 ) 70 Br 


70 Ge 


70 Br(p, 7 ) 71 Kr 


70 Br(p, 7 ) 71 Kr 




70 Se(p, 7 ) 71 Br 


70 Se(p, 7 ) 71 Br 




69 Se(p, 7 ) 70 Br 




71 As 


71 Br(p, 7 ) 72 Kr 


71 Br(p, 7 ) 72 Kr 


73 Se 


73 Kr(p, 7 ) 74 Rb 


73 Kr(p, 7 ) 74 Rb 


74 Se 


74 Kr(p, 7 ) 75 Rb 


74 Kr(p, 7 ) 75 Rb 
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Table 18 — Continued 



Lboi ope 


rusdCtlOIl ^illQlVlQUcll VclIlciClOnSJ 


riedction ^ivionxe odiioj 




74 Rb(p, 7 ) 75 Sr 

V-t 7 1 / 






73 Kr(p, 7 ) 74 Rb 




75 Br 


75 Rb(p, 7 ) 76 Sr 

\I 7 1 / 


75 Rb(p, 7 ) 76 Sr 

\Sr 7 1 / 


77 Kr 


77 Sr(p, 7 ) 78 Y 


77 Sr(p, 7 ) 78 Y 

\I 7 / / 




76 Rb(p, 7 ) 77 Sr 




78 Kr 


78 Sr(p, 7 ) 79 Y 


78 Sr(p, 7 ) 79 Y 

\-L 7 i / 




78 Y(p, 7 ) 79 Zr 




79 Kr 


79 Y(p, 7 ) 80 Zr 


79 Y(p, 7 ) 80 Zr 


81 Rb 


81 Zr(p, 7 ) 82 Nb 


81 Zr(p, 7 ) 82 Nb 

\I 7 1 / 




80 Y(p, 7 ) 81 Zr 




82 Sr 


82 Zr(p, 7 ) 83 Nb 


82 Zr(p, 7 ) 83 Nb 


83 Sr 


83 Nb(p, 7 ) 84 Mo 


83 Nb(p, 7 ) 84 Mo 




82 Zr(p, 7 ) 83 Nb 




84 Sr 


84 Nb(p, 7 ) 85 Mo 


84 Nb(p, 7 ) 85 Mo 


85y 


85 Mo(p, 7 ) 86 Tc 


85 Mo(p, 7 ) 86 Tc 


86 Zr 


86 Mo(p, 7 ) 87 Tc 


86 Mo(p, 7 ) 87 Tc 


87 Zr 


87 Tc(p, 7 ) 88 Ru 


87 Tc(p, 7 ) 88 Ru 




86 Mo(p, 7 ) 87 Tc 




88 Zr 


86 Mo(p, 7 ) 87 Tc 






87 Tc(p, 7 ) 88 Ru 




89 Nb 


86 Mo(p, 7 ) 87 Tc 






88 Tc(p, 7 ) 89 Ru 




90 Mo 


86 Mo(p, 7 ) 87 Tc 






85 Mo(p, 7 ) 86 Tc 




91 Nb 


85 Mo(p, 7 ) 86 Tc 
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Table 19. Summary of the most influential nuclear processes, as collected from Tables 
1-10. These reactions affect the yields of, at least, 3 isotopes when their nominal rates are 
varied by a factor of 10 up and/or down. See text for details. 



Reaction 


Models affected 


12 rV ^ „,\l6nia 

L>{a, 7) U 


rUo, 1YU4-B2, 1YU4-B4, iYU4-Bt> 


JNe(a, pj IN a 


T^n A Dib 
J\U4-B1 


25o:/„ „\28o 

oi(a, p) r 


K04-B5 


^Alya, pj di 


r Uo 


DyCt, pj L4 


rc; 


f{a, pj b 


rvU4-r>4 


30c/'„, ^33/^1 


T^n/l P/lb ~WC\A Rt;b 


Ul(_p, 7J Ar 


\CV\A P1 


32 c ^ _A36A V 

oyoz, 7 ) -".r 


T<Tl/1 RO 


im^q;, p j uu 


oUl , iYU4-rjD 


ou^p, 7 J Zjn 


r Uo 


59P,,/,, ^A60y„ 

uu(p, 7J Zn 


Qmb T^D/l R^ 
5U1 , rvU4-oO 


Ga(p, 7j be 


r08, K04-B1, K04-B2, K04-B5, K04-B6 


As^p, 7J oe 


x^ri/ib R1 T^ri/i Rob T^n/i RQb T^n/i R/i Tzn/i rc; izn/1 rk 
rvU4 , IYU4-B1, rvU4-BZ , rvU4-Bo , IYU4-B4, rvU4-B0, rvU4-B0 


Br^p, 7J Kr 


KU4-B i 


riD^p, 7J & r 


~Wf\A RO 


82 Zr(p, 7 ) 83 Nb 


K04-B6 


84 Zr(p, 7 ) 85 Nb 


K04-B2 


84 Nb(p, 7 ) 85 Mo 


K04-B6 


85 Mo(p, 7 ) 86 Tc 


F08 


86 Mo(p, 7 ) 87 Tc 


F08, K04-B6 


87 Mo(p, 7 ) 88 Tc 


K04-B6 


92 Ru(p, 7 ) 93 Rh 


K04-B2, K04-B6 


93 Rh(p, 7 ) 94 Pd 


K04-B2 


96 Ag(p, 7 ) 97 Cd 


K04, K04-B2, K04-B3, K04-B7 


102 In(p, 7 ) 103 Sn 


K04, K04-B3 


103 In(p, 7 ) 104 Sn 


K04-B3, K04-B7 


103 Sn(a, p) 106 Sb 


S01 b 
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a Reaction experimentally constrained to better than a factor of ~ 10 at XRB temperatures. 
See Section 5. 

b Reaction that affects the total energy generation rate by more than 5% at some time 
interval in this model, when its rate is varied by a factor of 10 up and/or down. See text 
and Table 20 for details. 
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Table 20. Nuclear processes affecting the total energy output by more than 5%, as well as 
the yield of at least one isotope, when their nominal rates are individually varied by a 
factor of 10 up and/or down, for the given model. See text for details. 



Reaction 


Models affected 


15 0(a, 7 ) 19 Ne a 


K04, K04-B1, K04-B6 


18 Ne(a, p) 21 Na a 


K04-B1, K04-B6 


22 Mg(a, p) 25 Al 


F08 


23 Al(p, 7 ) 24 Si 


K04-B1 


24 Mg(a, p) 27 Al a 


K04-B2 


26 3Al(p, 7 ) 27 Si a 


F08 


28 Si(a, p) 31 P a 


K04-B4 


30 S(a, p) 33 Cl 


K04-B4, K04-B5 


31 Cl(p, 7 ) 32 Ar 


K04-B3 


32 S(a, p) 35 Cl 


K04-B2 


35 Cl(p, 7 ) 36 Ar a 


K04-B2 


56 Ni(a, p) 59 Cu 


SOI 


59 Cu(p, 7 ) 60 Zn 


SOI 


65 As(p, 7 ) 66 Se 


K04, K04-B2, K04-B3 


69 Br(p, 7 ) 70 Kr 


SOI 


71 Br(p, 7 ) 72 Kr 


K04-B7 


103 Sn(«, p) 106 Sb 


SOI 



a Reaction experimentally constrained to better than a factor of ~10 at XRB temperatures. 
See Section 5. 
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Table 21. Comparison between theoretical rates from Tables 19 & 20, according to 
different rate compilations. The table shows, for a particular reaction, the maximum 
disagreement between our rate and the REACLIBvO rate over the temperature range 
covered in all our studies (~0.4 - 2.5 GK). All Q- values are given in keV. 



Reaction 


Q- value 


Q- value 


Q- value 


Maximum 




Audi et al. 2003a 


this work a 


REACLIBv0 a 


disagreement 


22 Mg(a, p) 25 Al 


3655(1) 


3655 [NS] 


3655 [NS] 


20% 


25 Si(a, p) 28 P 


6119(11) 


6122 [NS] 


6119 [NS] 


20% 


26 3Al(a, p) 29 Si 


4820.68(6) 


4823 [NS] 


4820 [NS] 


40% 


29 S(a, p) 32 Cl 


5306(50) 


5307 [NS] 


5306 [NS] 


x2 


30 S(a, p) 33 Cl 


2077(3) 


2076 [NS] 


2076 [NS] 


x2 


30 P(a, p) 33 S 


1521.36(34) 


1522 [NS] 


1521 [NS] 


20% 


31 Cl(p, 7 ) 32 Ar 


2422(50) 


2404 [Her95] 


2404 [Her95] 





32 S(a, 7 ) 36 Ar 


6640.76(14) 


6639 [NS] 


6641 [NS] 


x2 


56 Ni(a, p) 59 Cu 


-2411(11) 


-2413 [NS] 


-2411 [NS] 


x3 


57 Cu(p, 7 ) 58 Zn 


2277(52) 


2277 [ForOl] 


2277 [ForOl] 





59 Cu(p, 7 ) 60 Zn 


5120(11) 


5121 [NS] 


5120 [NS] 


10% 


61 Ga(p, 7 ) 62 Ge 


2442(149) b 


2427 [NS] 


2576 [Fis04] 


xlO 


65 As(p, 7 ) 66 Se 


2030(424) b 


2587 [NS] 


2433 [NS] 


10% 


69 Br(p, 7 ) 70 Kr 


2489(399) b 


2434 [NS] 


2579 [NS] 


10% 


71 Br(p, 7 ) 72 Kr 


4167(568) 


4525 [NS] 


4167 [NS] 


30% 


75 Rb(p, 7 ) 76 Sr 


4311 (38) 


4467 [NS] 


4311 [NS] 


10% 


Zr^p, 7 J 1Mb 


zUoo^oo ( ) 


zUbU [IN oj 


ZUOO [IN oj 


one/ 
20% 


84 Nb(p, 7 ) 85 Mo 


4513(409) b 


4510 [M] 


4513 [NS] 


40% 


84 Zr(p, 7 ) 85 Nb 


2946(297) b 


2109 [NS] 


2946 [NS] 


40% 


85 Mo(p, 7 ) 86 Tc 


1393(409) b 


1390 [M] 


1393 [NS] 


x4 


86 Mo(p, 7 ) 87 Tc 


1855(530) b 


1860 [M] 


1855 [NS] 


x3 


87 Mo(p, 7 ) 88 Tc 


2304(300) b 


2300 [M] 


2304 [NS] 


50% 


92 Ru(p, 7 ) 93 Rh 


2054(499) b 


2002 [NS] 


2054 [NS] 


20% 


93 Rh(p, 7 ) 94 Pd 


4467(566) b 


4365 [NS] 


4466 [NS] 


40% 


96 Ag(p, 7 ) 97 Cd 


3321(566) b 


3258 [NS] 


3321 [NS] 


x2 


102 In(p, 7 ) 103 Sn 


3554(318) b 


4095 [NS] 


3554 [NS] 


x2 


103 In(p, 7 ) 104 Sn 


4281(107) 


4241 [NS] 


4281 [NS] 


40% 


103 Sn(a, p) 106 Sb 


-5508(432) b 


-5508 [NS] 


-5508 [NS] 


10% 
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a Sources of the theoretical estimates: M=MOST Hauser-Feschbach code; NS=NON- 
SMOKER hauser-Feschbach code; For01=Forstner et al. (2001); Fis04=Fisker's shell model 
calculation -see REACLIBvO database; Her95=Herndl et al. (1995), shell model calculation. 

b Q- value estimated from systematics. See Audi et al. (2003a). 
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Fig. I. — Temperature versus time profiles corresponding to models K04 (Koike et al. 2004), 
SOI (Schatz et al. 2001), and F08 (Fisker et al. 2008). 
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Fig. 2. — Time evolution of the hydrogen mass fraction in models K04-B6 (low T) and K04- 
B7 (high T). Models are based on Koike et al. (2004) (see text for details). t=0 corresponds 
to the time at which T peak is reached. 
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Fig. 3. — Uncertainties in the final distribution of abundances resulting from the simultane- 
ous variation of all nuclear processes in bulk by factors ranging roughly between 0.1 and 10, 
for Model K04 (see Section 4.2 for details). 
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Fig. 4. — Same as Fig. 3, when all nuclear processes, except the 3a and all the /3-decay 
rates, are allowed to vary simultaneously by factors ranging roughly between 0.1 and 10, for 
Model K04 (see Section 4.2 for details). 
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Fig. 5.- 



Same as Fig. 3, for Model F08. 
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Fig. 6- 



Same as Fig. 4, for Model F08. 
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Fig. 7. — Correlation between the random enhancement factors applied to the 56 Ni(a, p) 
rate and the final 69 Ge normalized abundance, for Model K04. 
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Fig. 8. — Same as Fig. 7, for the strong correlation between the 69 Se(p, 7) rate and the final 
69 Ge abundance. 
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Fig. 9. — Same as Fig. 7, for the strong correlation between the 65 As(p, 7) rate and the final 
72 Se abundance. 



